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Carrier-free ^**^Po solution standards 
have been prepared and calibrated. The 
standards, which will be disseminated 
by the National Institute of Standards 
and Technology as Standard Reference 
Material SRM 4326, consist of (5.1597 
± 0.0024) g of a solution of polonium 
in nominal 2 mol • L"^ hydrochloric 
acid (having a solution density of 
(1.031±0.004) g • mL-i at 22 ''C) that 
is contained in 5 mL flame-sealed 
borosilicate glass ampoules, and are 
certified to contain a ^"''Po alpha-parti- 
cle emission rate concentration of 
(85.42 ±0.29) S-' • g-' (corresponding 
to a ^Po activity concentration of 
(85.83 ±0.30) Bq • g-^ as of the refer- 
ence time of 1200 EST 15 March 1994. 
The calibration was based on 4Tra liq- 
uid scintillation (LS) measurements 
with two different LS counting systems 
and under wide variations in measure- 
ment and sample conditions. Confirma- 
tory measurements by 2'ira gas-flow 



proportional counting were also per- 
formed. The only known radionuclidic 
impurity, based on a- and photon-emis- 
sion spectrometry, is a trace quantity of 
^^^^Po. The 2'***Po to 2(wpo impurity ratio 
as of the reference time was 0.00124 
±0.00020. All of the above cited uncer- 
tainty intervals correspond to a com- 
bined standard uncertainty multiplied 
by a coverage factor oi k = 2. Although 
^**^Po is nearly a pure a emitter with 
only a weak electron capture branch to 
^Bi, LS measurements of the ^*^Po a 
decay are confounded by an a transi- 
tion to a 2.3 keV (/'= 1/2") level in 
^^Pb which was previously unknown to 
be a delayed isomeric state. 
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1. Introductory Notes 

Solution standards of polonium isotopes have 
been and are important in a variety of applied mea- 
surement disciplines. At present, they are primarily 
of interest because of their use as calibration stan- 
dards for alpha-emission rate measurements, and 
as low-level tracers and separation yield monitors 
in radiochemical procedures that are employed 
with environmental and geophysical samples. 

For many years, preceding the past one or two 
decades, ^'^Po was unquestionably the polonium 
isotope most widely used in nuclear physics and ra- 
diochemistry research [1]. The reasons for this past 
popularity are evident: it could be prepared rela- 
tively easily and with virtually no radionuclidic im- 
purities by (d,n) reactions on natural monoisotopic 
^Bi, although historically it most often has been 



obtained by extraction from very aged natural sam- 
ples containing progenitors in the primordial ^^*U 
series; its half-life (138 d) is neither inconveniently 
too short nor too long for use in many physical and 
chemical studies; it has no intervening radioactive 
daughter products, decaying directly to stable ^^Pb; 
and its decay is by 100 % monoenergetic alpha 
emission. Over the years, however, its popularity 
has somewhat waned because of the advantages in- 
herent in having longer-lived standards for use 
within applied-measurement laboratories. 

More recently, in 1984, the National Bureau of 
Standards (NBS), now NIST, prepared, calibrated 
and disseminated carrier-free ^**^Po solution stan- 
dards. These longer-lived ^*Po standards (2.9 yr 
half-life) were issued as NIST Standard Reference 
Material SRM 4327. They consisted of approxi- 
mately 1.1 g of 1 mol • L"^ hydrochloric acid (HCl) 
contained in 2 mL flame-sealed borosilicate-glass 
ampoules, and were certified to contain a ^Po ac- 
tivity concentration of (76.6 ±1.1) Bq • g"^ on 20 
June, 1984 [2]. The above cited uncertainty interval 
corresponds to a combined standard uncertainty 
multiplied by a coverage factor of /: -3. The solu- 
tions contained a small quantity of contaminant 
^Po with a 2**^Po to ^Po impurity ratio of 0.0065 ± 
0.0026 at the reference time. The impurity ratio, of 
course, worsened with increasing age of the stan- 
dard. These SRMs had a great demand from labo- 
ratories throughout the world and the initial supply 
of several hundred sources was exhausted within a 
few years. 

Since the initial issuance of the SRM 4327 ^Po 
standards, considerable interest was expressed for 
having available similar standards of the even 
longer-lived 102 yr ^Po. In addition to the longer 
half-life, standards of the ^**^Po isotope also have 
the advantage of being capable of being prepared 
in a very pure state by ^Bi(p,n) reactions. Further- 
more, the impurity ratios for any slight traces of the 
shorter-lived ^Po or ^^^Po impurities will of neces- 
sity improve with time. 

The present work to prepare and calibrate a new 
series of ^Po solution standards follows a preced- 
ing companion study by Colle [3] in which he inves- 
tigated the long-term stability of very dilute, aged 
solutions of polonium. Prior to this earlier study, 
the integrity of polonium solutions stored over a 
decade or so was largely unknown. Considering the 
complex and often-perplexing chemistry of polo- 
nium [4-6], the previous study was considered to 
be advisable prior to issuing any further polonium 
standards, such as that for ^^o described herein. 
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2. Methods, Results, and Discussions 

2.1 Experimental Overview and Methodology 

2.1.1 Experimental Design The initial experi- 
mental design for the preparation and calibration 
of the standards was found to require considerable 
subsequent modification. The original design is 
outlined in the schema of Fig. 1. This design was 
followed, but, because of some difficulties and in- 
explicable early findings in the liquid scintillation 
(LS) measurements, a number of additional inde- 
pendent LS calibrations were also performed. 
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Fig. 1. Schema for preparation of the ^**^Po standards and 
counting sources used for the calibration. 

As indicated in Fig. 1, the ^^Po stock solution 
was transferred to a larger vessel (solution A) from 
which two LS samples were prepared. These solu- 
tion A samples were intended to be used to obtain 
an approximate assay of the ^^Po content prior to 
proceeding with the remaining scheme. Based on 
the findings of this assay, solution A could then be 
appropriately diluted to the starting point of solu- 



tion B. Four sets of counting sources that were to 
be used for the primary calibration were gravimet- 
rically prepared from solution B. The sources con- 
sisted of two sets of LS samples with two different 
quench conditions; three evaporated solid sources 
for use in a-emission impurity analyses and for 
confirmatory measurements by lira counting; and 
three evaporated solid sources for use in photon- 
emission impurity analyses. A master dispensing 
solution (M) was also prepared by a careful gravi- 
metric dilution of solution B. Approximately 200 
ampoules of the ^^Po standard were prepared from 
the master solution M. Four randomly selected am- 
poules were then used to prepare 4 sets of LS sam- 
ples to estimate between-ampoule variability. 

It was envisaged that the activity concentration 
of the ^^Po in the ampoule standards would be ob- 
tained from 4Tra LS measurements on the solution 
B samples after applying the known gravimetric di- 
lution factor to the assay results. Samples Bl 
through B3 were minimally quenched, and samples 
B4 through B6 were intended to be matched to the 
quench conditions of the M samples (see Fig. 1). 
Comparison of the LS results for the matched B 
and M samples would therefore verify the dilution 
factor. Confirmatory measurements were also per- 
formed by lua gas-flow proportional counting of 
the evaporated solid sources LI through L3. From 
the onset, it was considered that these solid source 
measurements would at best only serve as less ac- 
curate, but independent confirmations of the LS 
results. In contradistinction, the 1984 calibration of 
the previous ^^Po standards [2] was based on de- 
fined-geometry a-particle scintillation counting of 
evaporated sources, and used LS counting as an 
independent confirmatory measurement method. 
This difference in the respective calibration ap- 
proaches reflects, in part, the improvements in con- 
ventional LS standardization techniques over the 
past decade. 

2.1.2 Absolute Standardization by 4ira LS 
Counting and Problems Encountered with ^***Po 
Absolute standardizations of pure a-emitting ra- 
dionuclides by 4'n-a LS measurements are relatively 
straightforward.^ It was assumed, mistakenly, that 
this would be the case for this calibration. The LS 



^ Historical note: The use of LS counting as an "absolute" a 
standardization method was first proposed by Basson and Steyn 
[7] in 1954. Shortly thereafter, Basson [8] and Horrocks and 
Studier [9] demonstrated the feasibility and applicability of 4'n-a 
LS spectrometry for absolute standardization. Its use has since 
then become a fairly routine direct measurement methodology 
for a emitting radionuclides, including within the NIST Ra- 
dioactivity Group. 
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detection efficiency for a particles is virtually 100 % 
and independent of reasonable sample quenching 
conditions. Therefore, it was assumed that a cali- 
bration could be obtained by merely measuring and 
integrating the full-energy LS spectra for the ^^Po 
samples with appropriate corrections for back- 
ground subtractions, for extrapolations of the spec- 
tra to zero energy, for the presence of any 
impurities, and for decay of the ^^Po and any 
known impurities. The corrected LS count rate as a 
function of sample mass can also be extrapolated 
to zero sample mass (as a kind of quench correc- 
tion), but this is seldom necessary with a emitters 
since the sample mass dependency is usually found 
to be completely invariant of quench conditions. 
Lastly, it was recognized that the LS calibration 
would require a correction for the small (/ec = 
0.0048 ±0.0004) electron capture branch decay of 
2^^Po (refer to Fig. 2). 

The results of LS measurements based on the 
above described protocol evidenced a number of 
anomalies. 

First, all of the LS samples were measured with 
two different LS counting systems, and the two sets 
of corrected results exhibited substantial systematic 
differences. Depending on the specific sample in 
question, the differences in LS response between 
the two LS instruments ranged from about 0.5 % to 
several percent. This result would suggest a differ- 
ence in discrimination level or detection threshold 
for the two LS systems. Yet, this should not have 
been important for a nearly pure a emitter like 



^^Po. Such an instrument-dependent discrepancy 
had not been observed previously with other ra- 
dionuclides, including P" emitters like ^^Cl where 
the effect could be expected to be more pro- 
nounced [10]. 

Second, the magnitude of the differences be- 
tween the two instruments for a given sample ap- 
peared to be temporally dependent— i.e., 
dependent on the time of sample measurement 
This, of course, suggested the possibility of the de- 
tection efficiency (and or quenching) changing with 
time. Again, this should not have been important 
for ^^Po with an expected nearly quench-indepen- 
dent LS response. 

Third, early LS samples prepared from solutions 
A and B exhibited very rapid changes in instrument 
quench indicating parameters. This too should not 
have affected the LS a detection efficiency. 

Fourth, the dilution factors obtained from com- 
parison of the LS results with the solution B and M 
samples (using only the results from a single instru- 
ment) were not in as good agreement with the 
gravimetrically-determined dilution factor as that 
expected. The LS-determined dilution factors (for 
either instrument) differed from the gravimetri- 
cally-determined one by a few tenths of a percent. 
The observed range of differences was of course 
due to the temporal variations noted above. Typi- 
cally, for a single dilution, such dilution factor veri- 
fications differ by less than a few hundredths of a 
percent. 
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Fig. 2. Partial decay scheme for the ^**Po alpha and electron capture branch decays. 
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Fifth, and lastly, the LS spectra contained some 
structural details that were not understood or ex- 
pected. The magnitude of some of the observed 
differences and effects could not be accounted for 
by any known non-a decay mode in ^°^Po, such as 
that from electron-capture branching. 

Considering that these seeming measurement 
discrepancies were not experienced in many previ- 
ous LS standardizations by this laboratory, the 
causative culprit was initially suspected to be either 
some peculiarity of polonium chemistry, some type 
of LS cocktail instability, or some combination of 
the two. 

In an attempt to understand and reconcile the 
above observed discrepancies, a considerable num- 
ber of additional LS samples were prepared for 
further investigations. In large measure, this fur- 
ther work involved a systematic study of the LS re- 
sponses and spectral details as a function of 
different sample compositions and quenching con- 
ditions. The initial phases of the work was con- 
ducted in the period of March through June 1993 
followed by an eight-month abeyance. Further sys- 
tematic studies were continued and finally con- 
cluded in the period between February and April 
1994. 

2.1.3 ^^Po Decay Considerations The cause 
of the initially discrepant LS results was finally un- 
derstood during the 1994 period. It was deduced 
that they were not measurement discrepancies per 
se, but rather were a result of the initially chosen 
analysis method which was inappropriate when ap- 
plied to the measurement of ^**^Po. The problem 
arose from an inherent, but unknown aspect of 
^^Po decay: a sufficiently long, delayed isomeric 
state in ^^Pb. A partial decay scheme, based on the 
recent critical evaluation by Martin [11], is shown 
in Fig. 2. As indicated, the relative a branching to 
the two lowest energy levels in ^^'^Pb is largely un- 
known. Neither of the two most complete studies 
Qf 209pQ ^ decay [12,13] could resolve differences in 
an a transition to the ^^'^Pb ground state {P^SIT) 
or to the 2.3 keV (/^= 1/2") first excited state. The 
energy of the principal a branch (£« = 4.883 MeV) 
measured in these studies is assumed to corre- 
spond to the branch to the 2.3 keV level. The a 
transition to the ground state is thus 2 keV larger, 
i.e., £a- 4.885 MeV. The measured intensity for 
this combined a transition (/« = 0.9860 ±0.0008) is 
that for populating both levels. The assumed rela- 
tive populations of =»20 % to the ground state and 
~ 80 % to the 2.3 keV level are based on estimates 
by Martin [11] using the systematics of a-decay hin- 
derance factors. Irrespective of the absolute a tran- 



sition intensity to the 2.3 keV level, it seems clear 
that it is probably populated to a substantial extent. 
Furthermore, it is possible that the decay of this 
state is delayed and not in a prompt coincidence 
with the a feeding of it. The resulting l/2"->5/2" 
transition is also likely to be very highly converted 
with the emission of 2.3 keV conversion electrons 
(ce). If the lifetime of this 2.3 keV state is com- 
parable or longer lived than the pulse resolving 
time of the LS counting systems, than some frac- 
tion of these ce (and attendant 2.3 keV y rays, if 
any) would be detected independent of the feeding 
a-decay branch. The LS detection efficiency for the 
ce would also be dependent on sample quenching 
conditions and on the specific discrimination levels 
or detection thresholds for the LS instruments. 

Detailed examinations of the LS spectra as a 
function of quench conditions verified these con- 
jectures. The 2.3 keV ce transition was identified in 
spectra from both LS counters. A distinct low-en- 
ergy peak centered at about 3 keV was observed. 
Its large relative magnitude, ranging from just less 
than 1 % to nearly 6 % in various spectra, excluded 
the possibility of attributing the peak to other non- 
a decay modes, e.g., to Auger electrons from the 
weak electron capture decay branch. Its magnitude 
was also quench dependent, and was invariably 
larger in the instrument having a shorter coincident 
pulse resolving time. When this ce component was 
subtracted from the total LS response and spectra, 
one could obtain a virtually quench-independent a 
emission rate. Although the total LS detection effi- 
ciencies (which included the 2.3 keV component) 
for the two LS counting systems differed substan- 
tially, the subtracted a emission rates obtained 
from both instruments where in excellent agree- 
ment and invariant of quench conditions. 

Based on these findings, all of the original 1993 
results were also successfully re-analyzed and 
found to be in nearly perfect agreement with the 
1994 results. Hence, this standardization of ^^'^'Po 
was based on essentially two completely indepen- 
dent determinations that were separated by nearly 
a year in time, each of which was also based on the 
use of two different and independent LS counting 
systems. The standardization may be considered to 
be somewhat over-determined! It ultimately in- 
volved nearly 500 independent LS measurements 
with each instrument using over 40 different LS 
samples (covering a diverse range of sample com- 
positions and quenching conditions) which were 
obtained from 9 different ^Po solutions having 
gravimetrically related concentrations. 
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2.1.4 Measurement Uncertainty Treatments 

The specification of measurement uncertainties 
and the uncertainty analysis procedures used 
throughout this paper follow the normal conven- 
tions of the NIST Radioactivity Group. These con- 
ventions are wholly compatible with those adopted 
by the principal international metrology standard- 
ization bodies [14,15]. All individual uncertainty 
components are expressed in terms of estimated 
(experimental) standard deviations (or standard 
deviations of the mean where appropriate) or 
quantities assumed to correspond to standard devi- 
ations irrespective of the method used to evaluate 
their magnitude. All of these component quantities 
are designated as **standard uncertainties." A 
propagated uncertainty, termed a "combined stan- 
dard uncertainty," is expressed as what is assumed 
to be an estimated standard deviation which is 
equal to the positive square root of the total vari- 
ance obtained by summing all variance (square of 
the standard uncertainty) and covariance compo- 
nents, however evaluated, using the law of propa- 
gation of uncertainty for the specific mathematical 
function given by the model of the measurement 
procedure [14]. By recently established NIST pol- 
icy [15], the combined standard uncertainty is mul- 
tiplied by a "coverage factor" of A: = 2 to obtain an 
"expanded uncertainty" which is assumed to 
provide an uncertainty interval having a level of 
confidence of roughly 90 % to 95 %. For compara- 
tive purposes, it should be noted that previous 
SRM certificates issued by the NIST Radioactivity 
Group used comparably-based uncertainty cover- 
age factors of /: =3. This former practice was his- 
torically rooted and was assumed to provide 
certified uncertainty intervals with somewhat 
higher confidence levels, approaching 95 % to 
99%. 

Many of the component standard uncertainty es- 
timates, such as those for gravimetric determina- 
tions, detection efficiencies, and various correction 
factor and extrapolation procedures, are based on 
canonical values or on approximations by estab- 
lished procedures which were in turn determined 
by diverse experimental verifications and infer- 
ences and which were established by this labora- 
tory over many years of experience. 

2.2 Solution, Counting Source and Standard 
Preparations 

2.2.1 ^***Po Stock Solution The '"'Po stock so- 
lution used to prepare the standards has a some- 
what uncertain history. It was obtained by NIST 



from Oak Ridge National Laboratory (ORNL) in 
May 1989; and, based on partial records, appears 
to have been originally produced and radiochemi- 
cally separated by Mound Laboratory at least sev- 
eral years earlier [16]. The high degree of 
radionuclidic purity found in subsequent impurity 
evaluations suggests that it most probably was pre- 
pared by low-energy (p,n) reactions on natural 
monoisotopic ^**^Bi. The stock solution, contained 
and shipped in a 1 mL glass minivial, consisted of a 
nominal 100 kBq ^**^Po as polonium nitrate in about 
0.6 mL of 5 mol • L"' nitric acid (HNO3). 

2.2.2 Solutions A and B The stock solution 
was transferred and diluted with about 8 g of 
warmed 2 mol • L"^ HCl using an aspirating plastic 
pycnometer for the transfer and subsequent rinsing 
of the minivial. The resulting solution (labelled A 
in Fig. 1) was contained in a 20 mL wide-mouth 
glass vessel. Two gravimetrically determined ali- 
quants of the solution were used to prepare two LS 
samples (described below). Rudimentary LS assays 
of these samples indicated that most of the polo- 
nium was transferred; roughly 97 % of that said by 
ORNL to be contained in the minivial. If this had 
not been the case (as has sometimes been found in 
trying to remove polonium from glass vessels [3]), 
then more heroic efforts could easily have been 
made to recover the polonium remaining in the 
minivial and adding it to solution A. This fortu- 
nately was unnecessary in this particular transfer. 

Solution A was subsequently kept warmed and 
unstoppered over a period of about 60 h in an at- 
tempt to remove or decrease the HNO3 concentra- 
tion in the solution by evolving NO2. Warming of 
the solution was maintained with an infrared heat 
lamp at temperatures around 60 °C to 70 ''C. 
Higher temperatures were not employed in fear 
that the polonium might be hydrolyzed or 
volatilized and thereby lost from solution [4-6]. 
The resulting solution was labelled B. The activity 
concentrations of solutions A and B should hence 
be equal except for possible small evaporation 
losses during the warming. Solution B was then 
used to prepare the four sets of counting sources 
that are noted in Fig. 1 along with a master dis- 
pensing solution M. 

2.2.3 Gravimetric Determinations All of the 
solution aliquants used to prepare the counting 
sources and solution M were determined gravimet- 
rically with an estimated uncertainty, correspond- 
ing to an assumed standard deviation, of ± 10 ^g to 
±30 ^Jig. The mass determinations were made for 
the most part with single-pan suspension microbal- 
ances having internal balance reference weights. 
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For sample aliquants ranging from 50 mg to 100 mg 
for the counting sources, the relative standard un- 
certainties in the sample mass measurements thus 
are about ±0.01 % to ±0.03 % in general, and 
about ±0.06 % in worst cases, after appropriate 
gravimetric measurement air-buoyancy corrections, 
and considering the internal-balance-weight uncer- 
tainties, and typical mass measurement precision. 
Subsequent LS samples prepared from sealed am- 
poules of master solution M used sample masses 
typically ranging from 0.1 g to 1 g and thus had 
relative standard uncertainties of roughly ±0.05 % 
in worst cases. For solution M which had an ali- 
quant sample mass of greater than 7.5 g, the rela- 
tive standard uncertainty in the sample mass is less 
than 0.001 percent. In all cases, aliquant sample 
masses were redundantly and independently deter- 
mined from pycnometer-dispensed mass differ- 
ences and from contained masses obtained from 
the differences between the final sample masses 
and initial tare masses of the containing vessels. 
Differences in the two mass determinations are 
typically less than 0.05 percent for sample masses 
greater than 100 mg. This type of sample mass veri- 
fication was employed throughout the preparation 
of all samples given in this paper, and has mass 
uncertainties typical for the routine gravimetric 
procedures used in our laboratory. The uncertainty 
in the gravimetric dilution factor between solutions 
B and M is of course dominated by the uncertainty 
in the large HCl diluent mass which is in turn dom- 
inated by the uncertainties in the internal balance 
weights for the employed 2 kg single pan suspen- 
sion balance, and has an estimated overall relative 
standard uncertainty of about ±0.05 %. 

2.2.4 Master Solution M and Ampoules of the 
^***Po Standard In preparing solution M, solution 
B was diluted with 2 mol * L~* HCl by a gravimetri- 
cally determined factor of Z) = 147.243 ±0.074 
(about 7.5409 g of B in 1110.3 g of M). 

The solution M diluent was a carefully prepared 
HCl solution having a nominal acid concentration 
of 2.0 mol * L-' and density of (1.031 ±0.002) 
g • mL"* at 22 ''C. It was prepared from NIST high- 
purity, vacuum-distilled hydrochloric acid (nominal 
9.8 mol • L"^) in doubly-distilled water The uncer- 
tainty in the density is an estimated standard 
deviation based on independent gravimetric deter- 
minations of the total master solution mass and in- 
dividually dispensed masses in known dispensing 
volumes. This identical HCl solution was also used 
in the transfer of the stock ^*^Po solution to solu- 
tion A. It along with a blank doubly-distilled water 
solution was also used in preparing some of the 



subsequent LS samples. The blank water and HCl 
solutions were used to vary the HCl content in 
some of these LS counting sources, to vary or 
match sample compositions (and hence sample 
quenching) for the different LS sources, and to 
prepare matched blanks of nearly identical compo- 
sition for background subtractions. 

To prepare the ^Po standard sources, approxi- 
mately 5 mL of solution M was dispensed into each 
of 210 or so NIST 5 mL borosilicate-glass am- 
poules. Forty two, or about one out of every fifth 
ampoule that was filled, had been pre-weighed and 
then post-weighed to obtain an estimate of the dis- 
pensing precision and uncertainty in the contained 
solution mass. The ampoules were flame sealed af- 
ter filling and weighing. The mean mass of the dis- 
pensed solution in the ampoules was (5.1597 ± 
0.0012) g. The uncertainty given here is the stan- 
dard deviation of the 42 individually weighed 
masses of solution from which the mean was calcu- 
lated. It corresponds to a relative standard devia- 
tion of the mean of 0.0036 %. Figure 3 illustrates 
the dispensing precision as a function of ampoule 
filling. One may note that the mass variations ap- 
pear to be somewhat larger with the later fillings 
(higher ampoule numbers) compared to that at the 
beginning (lower ampoule numbers). It is quite 
probable that this is not actually due to poorer dis- 
pensing precision into later ampoules, but rather is 
due to greater imprecision in the mass measure- 
ments at later times. One can speculate that this 
arose as the result of increasing imprecision in the 
balance readings introduced by fatigue of a labora- 
tory worker who had to solely fill and seal 200^ 
ampoules, along with pre- and post-weighing a sub- 
stantial fraction of them, all of which was per- 
formed over the course of 5 h. 
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Fig. 3. Dispensed solution masses of the ^*^Po standards as a 
function of ampoule number filling. The solid line is the mean 
of 42 determinations. The broken lines are upper and lower 
uncertainty limits for plus and minus two standard deviations 
about the mean. 
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2.2.5 2ira- and Photon-Emission Point Sources 

Three 2'Tra point sources (labelled LI through L3 
in Fig. 1) were prepared by depositing gravimetri- 
cally determined aliquants of solution B onto 2.5 
cm diameter polished stainless steel (SS) disks. 
The deposits ranged from 39 mg to 51 mg of solu- 
tion B, were dispersed and seeded with Ludox (E.L 
DuPont de Nemours, Wilmington, DE),^ and then 
were evaporated to dryness without heating in am- 
bient air. The acid content of solution B, however, 
seriously corroded and pitted the surfaces of the SS 
source mounts which was probably intensified by 
remaining traces of HNO3 in the solution. The 
sources were found to be unusable for quantitative 
a counting or a spectroscopy. It is also suspected 
that polonium may have been partially lost from 
the sources since the very low observed a emission 
rates (27 % to 36 % of that expected) seems too 
extreme to be accounted for solely by source self- 
absorption losses. Two additional sources (not 
shown in Fig. 1) were prepared at a later time and 
labelled LXl and LX2. They consisted of 207 mg 
and 66 mg, respectively, of gravimetrically deter- 
mined aliquants of solution M (obtained from a 
randomly selected ampoule of the standard solu- 
tion) which were deposited onto 2.5 cm diameter 
platinum disks. The somewhat more massive solu- 
tion deposits were used because of the lower activ- 
ity concentration of solution M. The deposits 
covered a surface area of roughly 4.5 cm^ and 1.5 
cm^ and were similarly evaporated to dryness in air. 
Sources LXl and LX2 were used both for confir- 
matory measurements by gas-flow proportional 
counting and for spectroscopic analyses of a-emit- 
ting impurities. 

Three point sources (labelled Gl through G3 in 
Fig. 1) for y- and x-ray measurements were also 
prepared. They consisted of 32 mg to 82 mg of so- 
lution B deposited onto annular source mounts 
backed with «7 mg • cm"^ plastic tape (a glued 
polyester film). After the deposits air dried, the 
sources were covered and sealed with an identical 
layer of the tape. 

2.2.6 LS Counting Sources The sample com- 
positions for the LS counting sources are summa- 
rized in Table 1. Excepting a few samples that were 
prepared for special spectral investigations (and 
which will be noted later), the LS counting sources, 



^ Certain commercial equipment, instruments, or materials are 
identified in this paper to foster understanding. Such identifica- 
tion does not imply recommendation or endorsement by the 
National Institute of Standards and Technology, nor does it im- 
ply that the materials or equipment identified are necessarily 
the best available for the purpose. 



in all cases, consisted of an appropriate aliquant of 
one of the ^^Po solutions and approximately 10 mL 
or 14 mL of a prepared commercially-available LS 
cocktail mix contained within nominal 22 mL LS 
counting vials. The cocktail was "Ready Safe" 
(Beckman Instruments, Fullerton, CA), a prepared 
scintillation fluor in a polyarylalkane (i.e., an alky- 
lated biphenyl) solvent having surfactant agents 
[17]. Some of the sources also contained respective 
proportions of a blank doubly-distilled water solu- 
tion and the carefully prepared blank HCl solution. 
The water and HCl solution additions were done, 
as noted earlier, to vary the LS sample composi- 
tions to study sample quenching effects. In all 
cases, the sources were compared with matched 
blanks of nearly identical composition for back- 
ground subtractions. The 1993 LS samples (Al 
through N6 in Table 1) were prepared in high-den- 
sity polyethylene LS vials, whereas the 1994 
sources were contained in conventional glass LS 
vials. The change in vials was made to exclude the 
possibility that the initially inexplicable LS results 
were due to the use of plastic vials. 

Table 1 summarizes the specific ^^Po solution 
used to prepare the sources; the range for the mass 
rris of the ^^Po solution aliquant used in the source; 
the range for the total source mass nitot (which is 
the sum of nis and the masses of the scintillation 
cocktail, nil, and any blank water, m*, or blank 
HCl, wia, used in the source); the range for the 
aqueous phase content in the source, expressed as 
a percentage (paq = 100 • [ms4-mw + Wa]/mtot); and 
the range for the nominal normality of the acid 
content in the aqueous phase (^aq = 2.0 [ms + ma]/ 
[ms + mw+w2a]). Although the values given in Table 
1 are summarized in terms of approximate ranges, 
it should be recalled that each of the ms, Ww, Wa, 
and nit masses in every sample was determined 
with an estimated relative standard uncertainty of 
±0.05 %. 

Samples Al and A2 and Bl through B3 were 
prepared nearly identically with about ms=52 mg 
to 65 mg of solution A or B in about mr=^10 g of 
cocktail. Corresponding matched blanks contained 
ma— 51 mg to 60 mg of the blank HCl solution in 
mr^^^lO g of cocktail. Samples B4 through B6 were 
prepared with W2s=50 mg and W2w=l g in W2r=10 g. 
The additional water in these samples would essen- 
tially match the total sample masses (or volumes) 
of the MAI through MD3 sample series. Use of the 
blank water also more closely matched the quench- 
ing conditions in the two sample sets. Following an 
initial series of LS measurements on samples B4- 
B6, the samples were adjusted to the approximate 
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Table 1. LS 


sample compositions 










Sample 


209p^, 


ms Range*' 


mtot Range*" 


/?aq Range"* 


A^aq Range*^ 


Identities 


Solution* 


(g) 


(g) 


(%) 


(mol-L-^) 


A1-A2 


A 


0.05-0.07 


10.10-10.08 


0.54-0.65 


2.00 


B1-B3 


B 


0.05-0.06 


10.09-10.10 


0.52-0.55 


2.00 


B4-B6 


B 


0.05 


10.09-11.14 


9.44-9.84 


0.097-0.101 


B4x-B6x 


B 


0.05 


11.35-11.42 


11.5-12.1 


1.98-2.08 


MA1-MA3 


M (amp. 39) 


1.01-1.04 


11.05-11.08 


9.11-9.34 


2.00 


MB1-MB3 


M (amp. 93) 


0.99-1.04 


11.03-11.08 


8.97-9.36 


2.00 


MC1-MC3 


M (amp. 129) 


0.99-1.02 


11.03-11.06 


8.97-9.23 


2.00 


MD1-MD3 


M (amp. 205) 


1.00-1.01 


11.04-11.05 


9.08-9.15 


2.00 


N1-N6 


M (amp. 47) 


0.12-1.01 


11.11-11.13 


9.03-9.21 


0.225-2.00 


P1-P5 


M (amp. 40) 


0.11-0.97 


10.98-11.03 


8.87-9.15 


0.209-2.00 


01-04 


M (amp. 29) 


0.05-0.71 


10.06-11.03 


0.534-9.24 


0.225-2.00 


05-08 


M (amp. 29) 


0.13-1.00 


14.11-14.99 


0.892-6.70 


2.00 



" See Fig. 1. 

** Approximate mass of ^**^Po solution used in preparing the sample. 

*" Approximate total mass of sample composed of the sum of m^, Wr, Ww, and m^. Refer to text. 

** Approximate water content of the sample on a mass percentage basis. 

' Approximate acid content (normality) in the sample aqueous phase. 



A^aq=2 mol • L'^ acid normality of the M samples. 
This was achieved by the addition of approximately 
0.25 g to 0.3 g of the very pure NIST concentrated 
HCl. These adjusted samples are identified in 
Table 1 as samples B4x through B6x. The corre- 
sponding matched blanks were similarly adjusted. 
Four randomly selected sealed ampoules of the 
^°^Po solution standard were used to prepare the 12 
M series samples shown in Fig. 1. Each contained 
ffts^'l g in Wr=10 g. All further samples (not shown 
in Fig. 1) were also prepared using the sealed am- 
poules. Sample series Nl through N6, PI through 
P5, Ql through Q4, and Q5 through Q8 (along 
with all of their corresponding blanks) were pre- 
pared with varying quantities of Ws, Ww, Wa, and Wr 
to achieve varying ;?aq and A^aq sample compositions 
that could be used to systematically study cocktail 
stability, sample quenching, and spectral effects. 
Compositions of samples used in special spectral 
investigations are given later. 

2.3 LS Measurement Systems 

2.3.1 Instruments The two LS counting sys- 
tems were: (1) a Beckman LS7800 model LS coun- 
ter equipped with two Hamamasutu R331-05 
photomultiplier tubes (PMT) operating in a coinci- 
dence mode, a logarithmic pulse amplifier coupled 
to an analog-to-digital converter (ADC) for spec- 
tral pulse-height analysis, and an external *^^Cs 
source for Compton-edge (Horrocks number) 
quench monitoring [18]; and (2) a Packard Tri- 
Carb 2500TR LS analyzer employing two matched 



high performance PMT also operating in coinci- 
dence, but with linear amplification for ADC 
pulse-height spectral analysis, and an external low- 
energy y-ray ^^^Ba source that is used to obtain an 
instrument-provided "transformed Spectral Index 
of the External Standard (tSIEy quench indicating 
parameter [19]. 

Resolving times for the two instruments are 
slightly different [20,21]. Their relative magnitudes 
were particularly of interest in terms of interpret- 
ing the data in regard to the existence of the de- 
layed isomeric state in ^^^Pb. The resolving times of 
the coincidence gates for simultaneous detection of 
singles PMT events were comparable: Tc = 22 ns 
and Tc= 18 ns in the Beckman and Packard instru- 
ments, respectively. Of greater interest, are the re- 
solving times (i.e., equivalent deadtimes) for 
detection of coincident pulses. These were: t=5 |jls 
to 33 |JLS in the Beckman, and t=12 |jls in the Pack- 
ard. The latter instrument has a fixed (i.e., "non- 
extending") resolving deadtime. The r ranges for 
the former result from the resolving time's depen- 
dence on pulse height. The spectrometer would of 
course spend the vast majority of its time counting 
the principal ^^Po a. Using the pulse heights for 
the dominant a peak (i.e., its channel number loca- 
tions), the typical resolving times for the Beckman 
instrument were r = (15 ± 4) |jls. 

2.3.2 Quench Indicating Parameters The 
Horrocks number (H^) quench indicating parame- 
ter used with the Beckman system is based on the 
downward spectrum shift of the Compton edge of 
the external *^^Cs y-ray standard with increasing 
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sample quenching. The parameter corresponds to 
the spectral channel number shift between the 
quenched sample and an unquenched blank sam- 
ple. The channel number shift (c2— ci) is, in this 
case, equal to a logarithmic energy ratio \og(E2/Ei). 

The Packard system's quench indicating parame- 
ter tSIE is based on a proprietary mathematical 
transform of the energy distribution of the "^Ba 
generated Compton spectrum. The transform is 
said to be used to correct for spectral distortions 
arising, for example, from wall effects, volume vari- 
ations, and color quenching. The parameter con- 
sists of a relative quenching scale in which 
unquenched samples correspond to tSIE = 1000. 

The quench parameters H"^ and tSIE are in- 
versely correlated and scale proportionally with 
variations in the LS sample composition that effect 
the chemical quenching. Figure 4, for example, il- 
lustrates the correlation between the two parame- 
ters for the Nl through N6 samples prepared and 
measured in 1993. The slope d(r5/£)/d(//'*') corre- 
sponds to a decrease of about 3.5 units of tSIE for 
every unit ofH^. A typical and representative illus- 
tration of the effect of LS sample composition 
changes is shown in Fig. 5. It shows the relative 
change in the two quench parameters as a function 
of the acid content A^aq in the PI through P5 sam- 
ples which were prepared and measured in 1994. 
The values given in Fig. 5 are ratios of //'*'/ 
(93.65 ±0.69) and (439.S± 12)/tSIE that consist of 
normalizations to the quench parameters for the 
least quenched P5 sample. The quench parameters 
similarly scale with other LS sample composition 
variables such as nis or/?aq (see Table 1). 
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least quenched sample. Refer to text for normalizations. 



2.4 1993 4ira LS Measurement Results 

2,4,1 Analysis Procedures The results given 
here are based on the 1994 re-evaluation of the 
originally collected 1993 data. This re-evaluation 
took account of the confounding by the 2.3 keV 
delayed isomeric state in ^"^Pb. Two corrected LS 
counting rate concentrations were evaluated. The 
first Rac is that obtained from integrating the full- 
energy LS spectra, while the second Ra is that ob- 
tained after subtraction of the component 
responses from the 2.3 keV ^"^Pb transition. The 
integrated LS spectra energy regions, in terms of 
approximate beta energy, for the two are roughly 
from about AE^c = keV to 2000 keV for R^c, and 
from about AEa^30 keV to 2000 keV for /?«. The 
two counting rate concentrations are given by 



Rac — 



yL^ae i-^baej 



(e-'2^^+/e"'2(«^rm 



Ra = 



(Cg — Cba) 



(e-^^'^^+/e"'^*«0^m 



where Cae = gross integral counts in the region AE^c 
obtained from a spectrum of a given 
LS source at measurement time 7; 
Cbae = the corresponding gross integral counts 
in region AE„c for background ob- 
tained from a spectrum of a matched 
blank LS sample; 



10 
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Ca = gross integral counts in the region AEa 

of the source spectrum; 
Cba = the corresponding gross integral counts 
in the region AEa for the background 
spectrum; 
t = counting (live) time interval for mea- 
surement of both the source and 
matched blank; 
m =mass of the ^^Fo solution used in the 

LS source; 
A209 = ^^Po decay constant; 
A208 = ^"^Po decay constant; 
T = the decay time interval from the mid- 
point of the measurement time interval 
to a common reference time; 
I _208pQ ^Q 209pQ jjjipurity ratio in terms of 
their respective a emission rates at the 
reference time. 
The difference Re = (Rac-Ra) is, of course the ap- 
proximate LS counting rate concentration contri- 
bution due to the ^°^Pb 2.3 keV transition 
(presumed to be due predominantly to ce). 

Two very small corrections, both less than 0.1 % 
and of opposing magnitude, were later applied to 
Ra to obtain the ^**^Po solution standard a-particle 
emission rate concentration ^a. The first correction 
factor kzE consisted of an extrapolation of Ra to 
zero energy by extrapolating Ra under the Re por- 
tion of the spectra. The second correction factor 
kEc was to remove LS response contributions aris- 
ing from the EC branch of ^Po decay from Ra. 

The lower limits of the energy regions AE^ used 
to evaluate the 1993 data were slightly different for 
the two LS systems: for the Beckman instrument, 
AEa^30 keV (A and B sample series) and AEa ^45 
keV (M and N sample series); for the Packard, 
AEa^l9 keV (all sample series). Evaluation of the 
later 1994 LS data also used slightly different inte- 
grating regions (see Sec. 2.5). The conversions of 
spectral channel numbers to a beta energy scale 
are, at best, only very approximate. First, one must 
recognize that the light energy transfer mecha- 
nisms in the scintillator/solvent mixtures are very 
non-linear at low energies. Second, the channel 
number to energy conversions are dependent on 
sample quenching. Third, and lastly, it is difficult to 
directly compare the spectra from the two LS sys- 
tems since the Beckman instrument with its loga- 
rithmic amplification gives ADC channel numbers 
corresponding to a logarithmic energy scale while 
the conversion of channel number to energy is 
nominally linear with the Packard. 

A comparison of the ^**^Po LS spectra obtained 
with the two instruments is shown in Fig. 6. The 
spectra were obtained in 1994 with LS sample P3. 



The LS response from the dominant 4.88 MeV a of 
^^Fo appears as a peak at around roughly 300 to 
500 keV in the corresponding beta energy scale. 
This results from a particles (with their relatively 
short ranges) having scintillation yields that are 
typically a factor of 8 to 12 lower than that from 
beta emitters. The LS response from the ^°^Pb con- 
version electron, centered at approximately 3 keV, 
is readily apparent in both spectra. The relative 
magnitudes of the conversion electron (Re) and al- 
pha peaks (Ra) shown in the spectra may deceive 
on first appearance. One should note that the ordi- 
nate axis is in units of counts per second per keV, 
and that the abscissa is a logarithmic energy scale. 
The ratios of RJRa in the two spectra of Fig. 6 are 
actually 0.026 and 0.044 for the Beckman and 
Packard instruments, respectively. The large differ- 
ence in the RJRa values may be attributed to differ- 
ences in detection and resolution of the ^^Pb 
conversion electron (ce) by the two instruments. 
One may also note that in the region below 10 keV, 
the spectrum obtained with the Packard instrument 
has only 9 or 10 channels of data while that from 
the Beckman instrument has over 50. Values of Ra 
for the two spectra shown in Fig. 6 are in better 
than 0.2 % agreement. In the region above about 
550 keV, both spectra have net counting rates after 
background subtraction of virtually zero. The spec- 
tra in the region from about 30 keV to 100 keV are 
very flat and have net counting rates that are also 
very small which permits a relatively easy extrapo- 
lation ofRcc below the Re peak. Additional spectral 
analysis considerations are treated below (Sec. 2.6). 
The 1993 LS measurements were based on aver- 
ages of 5 to 10 replicate measurements of each 
sample interspersed with background measure- 
ments of blank samples. Each measurement ranged 
from 30 min to 120 min in duration. Because the 
initial results were not understood many of the 
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samples were remeasured from 2 to 5 more times 
(with each subsequent sequence consisting of 5 to 
10 measurements) over periods of nearly 100 d. 
The averaged results for R^c and R^ were decay 
corrected to the common reference time of 1200 
EST 15 March 1994 for direct comparison with the 
1994 LS results. Values of T used for the decay 
corrections therefore often exceeded a year. The 
value of 7=0.00124 used for the impurity correc- 
tion was independently obtained from a spec- 
trometry (see Sec. 2.8). 

2.4.2 A and B Series Samples All of the A 
and B series LS samples exhibited large and not 
very systematic variations in both quench indicating 
parameters and in values of /?ac. The measurement 
results for sample Al as a function of measurement 
time, shown in Fig. 7, are typical. Each plotted 
value in Fig. 7 is based on the mean of 5 to 10 
measurements with the Beckman instrument. Over 
a period of approximately 100 d, the quench 
parameter //'*' varied from about 75 to over 160, 
with the largest rate of change occurring in the first 
10 or so days. The Packard's corresponding tSIE 
for sample Al, over a shorter 55 d period, ranged 
from 386 to 261. For 5 to 10 replicate measure- 
ments, these quench parameters //* and tSIE were 
typically determined with a precision having a rela- 
tive standard deviation of the mean (5m) of 0.1 % to 
0.3 %. Along with the change in //*, the mean val- 
ues of Rctc over the same time —100 d interval de- 
creased by over 2 %. The counting results for every 
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centrations /?„ obtained from integrating the full-energy LS 
spectra; the lower curve is for /?« and is that obtained after 
subtraction of the component LS responses from the 2.3 keV 
2**5pb transition. 



A and B series sample measured with either instru- 
ment exhibited similar changes but of slightly dif- 
fering magnitudes. 

Later LS samples prepared from the standard 
solution ampoules did not undergo such large tem- 
poral variations in quench parameters and Km- It 
may have occurred with only the A and B series 
samples because of larger quantities of unexpelled 
HNO3 in these more concentrated solutions. The 
HNO3 content of the latter master solution used to 
prepare the ampoules was diluted by a factor of 
nearly 150. Nitric acid could very well have a more 
deleterious effect on the LS cocktails than the HCl 
diluting solution and could therefore possibly ac- 
count for the differences observed with the various 
LS samples. Solution A had a HNO3 normality of 
about 0.35 mol • L~\ and solution B may have been 
comparable if little of HNO3 was effectively re- 
moved by the gentle heat treatment. 

Although Rm was found to vary with time, with 
sample quenching conditions, and with the instru- 
ment used to perform the measurement, the values 
of Ra were invariant. Figure 7 illustrates this nicely. 
The mean values of Ra for sample Al over the 
same time interval are constant and independent of 
the sample quenching changes. This demonstration 
is achieved even more powerfully in Figs. 8 and 9 
which aggregate all of the measurement results on 
the series A and B samples. Each plotted value in 
the figures represents the mean of 5 to 10 measure- 
ments of either Rac (closed circles) or Ra (closed 
triangles). The multiple values shown for a given 
sample are the mean values for additional 5 to 10 
measurement sequences obtained at later times 
(and with changing quench parameters). The illus- 
trated error bars correspond to the calculated 
standard deviations for the 5 to 10 replicate mea- 
surements. They correspond to 5m values ranging 
from 0.02 % to 0.06 %. Figures 8 and 9 give the 
results obtained with the Beckman and Packard in- 
struments, respectively. Examinations of the fig- 
ures indicate that although the values of Rac vary 
substantially between instruments, between sam- 
ples, and at different times for any given sample, 
all of the Ra values are essentially constant. This is 
exactly what is to be expected: invariant and 100 % 
a detection efficiencies resulting in Ra values that 
are independent of quenching conditions or 
changes; very quench dependent and instrument 
dependent Rt values for the low energy ^"^Pb ce; 
and a greater R^ response in the instrument with 
shorter pulse resolving time. 

The mean values of Ra across all samples and 
measurements are 12560 s"^ • g"^ for the Beckman 
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obtained with the Packard LS system for the A and B series samples. Refer to Fig. 8 
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(Fig. 8) results, and 12553 s"^ • g"^ for the Packard 
(Fig. 9).^ The two respective mean values have s^ 
values of 0.015 % and 0.017 % based on totals of 
197 and 149 replicate measurements. The mean 
value for all 346 measurements obtained with both 
instruments is 12557 s"^ • g~^ with 5m = 0.016 %. On 
application of the gravimetrically determined dilu- 



^ Because all of the counting rate concentrations that follow, 
such as for /?„, R^^, and /?e, include the concept of "counts" in 
their definitions, their proper units are s"^ • g"* rather than the 
traditionally used "counts per second per gram." 



tion factor D= 147.243, the corresponding count- 
ing rate concentration for the ^^Po solution stan- 
dard is RJD =85.28 s"^ • g-\ 

2,4,3 M Series Samples The differences in 
observed R^ values obtained with the two LS sys- 
tems, due to differences in the detection and reso- 
lution of the ^^Pb ce, is even more clearly 
demonstrated in the results for the 12 M series of 
samples. The samples (three prepared from each 
of four randomly selected sealed ampoules) were 
of closely matched composition and had average 
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quench parameters ranging from H^ = 134 to 135 
and tSIE =2S7 to 291 which did not appreciably 
alter over a two week period. The observed values 
for Rae and Ra obtained with the two LS counters 
are illustrated in Fig. 10. Each plotted determina- 
tion is based on 10 replicate measurements of each 
sample on each instrument. The grand means Km 
(Beckman) = (86.89 ±0.26) s"^ - g"^ and Rac (Pack- 
ard) = (87.75 ± 0.24) S-' • g-^ differ by about 1 %, 
while the two grand means Ra (Beckman) = 
(85.43 ± 0.25) s"^ • g"^ and Ra (Packard) = (85,43 ± 
0.23) s"^ • g"* are in agreement to four significant 
figures. The above uncertainty intervals given for 
these means are standard deviations of the means 
(sm) with 119 degrees of freedom. These results 
also demonstrate the excellent between-sample 
and between- ampoule precision. Table 2 (for the 
Beckman data) and Table 3 (for the Packard) sum- 
marizes the precision estimators 5m for 10 replicate 
measurements on any one sample, and for the 3 
means of samples prepared from any one ampoule. 
All of the tabulated mean Ra values are statistically 
equivalent, independent of the instrument used, 
specific sample, or ampoule used to prepare the 
samples. 



2.4.4 N Series Samples The systematic varia- 
tion of Re (obtained from the differences Rac—Ra) 
with sample quenching is exhibited in the results of 
Figs. 11 and 12 for the N series of samples. These 
data are based on 5 replicate measurements of 
each sample with the Beckman instrument, and 10 
measurements on each sample with the Packard. 
Their quench parameters were previously provided 
in Fig. 4. Sample quenching was systematically 
varied by changing the sample compositions, in this 
case by varying the mass of ^^Po solution m^ used 
to prepare the samples. Total sample masses and 
aqueous content p^q were kept nearly constant by 
the additions of blank water (see Table 1). The ex- 
trapolations to zero solution mass (ms = 0) of Figs. 
11 and 12 could equally have been performed using 
the acid content A^aq variable. Figures 11 and 12 
also illustrate the replicate measurement precision 
for these LS measurements since all measurement 
results on each sample are plotted along with their 
respective error bars that correspond to the total 
statistical '^counting errors.'"* The precision is typi- 
cal of that obtained with all samples prepared from 
the sealed ampoules. 
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Fig. 10. Mean LS counting rate concentrations R^c and /?« obtained with the two LS systems for the 
M series samples. Closed squares are for R^ (Packard); closed circles for R^ (Beckman); open 
squares for R^ (Packard); and open circles for Ra (Beckman). 



* The total "counting error" devolves from the assumption that the 
radioactive decay process is Poisson distributed and that its variance 
(N) is equal to the mean number of total observed counts (N). The 
propagated Sj, is expresse d as a relative standard deviation and is 
merely given by VA^s + ^(A^s - A^b) where A^s is the total number of 
counts obtained in measuring the sample and A^b is the respective 
total number of background counts, and where A^s and A^b are ob- 
tained over counting time intervals of equal duration. The standard 
deviation in the counting rate concentration R (due to the counting 
errors) is then the product Su(R) =sji. 
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Table 2. LS measurement results for the M sample series obtained with the Beckman 
counting system 



Sample 


Number of /i " 


Ra mean** 


s"" 


5.^ 


identity 


samples 


(s-'-g-O 


(s-^*g-') 


(%) 




Replicate measurements on single 


samples 




MAI 


1 10 


85.414 


0.277 


0.103 


MA2 


1 10 


85.452 


0.222 


0.082 


MA3 


1 10 


85.617 


0.275 


0.102 


MBl 


1 10 


85.554 


0.347 


0.128 


MB2 


1 10 


85.437 


0.243 


0.090 


MB3 


1 10 


85.353 


0.125 


0.046 


MCI 


1 10 


85.370 


0.264 


0.098 


MC2 


1 10 


85.348 


0.251 


0.093 


MC3 


1 10 


85.282 


0.169 


0.063 


MDl 


1 10 


85.431 


0.164 


0.061 


MD2 


1 10 


85.416 


0.258 


0.096 


MD3 


1 10 


85.545 


0.224 


0.083 




Replicate measurements on 


samples from 


a single ampoule 




MA1-MA3 


3 30 


85.494 


0.266 


0.057 


MB1-MB3 


3 30 


85.448 


0.260 


0.056 


MC1-MC3 


3 30 


85.333 


0.227 


0.049 


MD1-MD3 


3 30 


85.464 


0.219 


0.047 


All measurements 


MA1-MD3 


12 120 


85.435 


0.253 


0.027 



" Total number of measurements on all samples in the series. 

^ Mean of the ^**^Po a counting rate concentration, with ^"**Po impurity and decay correc- 
tions to the reference time 1200 EST 15 March 1994 (see Sec. 2.4.1). 
*^ Standard deviation for the n determinations. 
^ Relative standard deviation of the mean in percent for the n determinations. 
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Fig. 11. LS counting rate concentrations Rac (closed triangles) and Ra (closed circles) obtained with 
the Beckman instrument for the N series samples as a function of m, (and sample quenching). The 
solid lines are linear regressions fitted to the data. 
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Table 3. LS measurement results for the M sample series obtained with the Packard 
counting system 



Sample 


Number of n * 


Ra mean** 


5*= 


5™** 


identity 


samples 


(s-^-g-^) 


(s-^-g-0 


(%) 


Replicate measurements on single samples 


MAI 


1 10 


85.563 


0.197 


0.073 


MA2 


1 10 


85.370 


0.190 


0.070 


MAS 


1 10 


85.321 


0.139 


0.052 


MBl 


1 10 


85.446 


0.241 


0,089 


MB2 


1 10 


85.459 


0.260 


0.096 


MB3 


1 10 


85,400 


0.229 


0,085 


MCI 


1 10 


85.304 


0.200 


0.074 


MC2 


1 10 


85.447 


0.317 


0.117 


MC3 


1 10 


85.562 


0.275 


0.102 


MDl 


1 10 


85.394 


0.215 


0.080 


MD2 


1 10 


85.445 


0.202 


0.075 


MD3 


1 10 


85.478 


0.233 


0.086 




Replicate measurements on 


samples from a 


single ampoule 




MA1-MA3 


3 30 


85.418 


0.201 


0.043 


MB1-MB3 


3 30 


85.435 


0.237 


0.051 


MC1-MC3 


3 30 


85.438 


0.280 


0.043 


MD1-MD3 


3 30 


85.439 


0.212 


0.045 


All measurements 


MA1-MD3 


12 120 


85.433 


0.232 


0.025 



■ Total number of measurements on all samples in the series. 

^ Mean of the ^'*'Po a counting rate concentration, with ^****Po impurity and decay correc- 
tions to the reference time 1200 EST 15 March 1994 (see Sec. 2.4.1). 
^ Standard deviation for the n determinations. 
'^ Relative standard deviation of the mean in percent for the n determinations. 
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The weighted ;^^-minimized (least-squares) re- 
gression lines shown in Fig. 11 (for the Beckman 
results) have fitted parameters of: 

/?ac = (86.769±0.083)s-^-g-i 

-[(0.55±0.14)s-^-g-^]Ws, 

and 

/?a=(85.451±0.083)s-i-g-i 

+ [(0.077±0.14)s-^-g-2]ms. 

The linear fits to the data were performed with 
weighting factors equal to the reciprocal variances 
ysp^(R). The fits were performed using all 30 mea- 
surement values. The standard deviations on the 
fitted parameters thus have v = (30-2) degrees of 
freedom. The magnitude and uncertainty on the 
slope d/?a/dms = (0.08 ±0.14) s"* • g"' clearly indi- 
cates that Ra is invariant of m& and of accompany- 
ing quench changes. Similarly, the Packard results 
of Fig. 12 are: 

/?„e = (88.807±0.097)s-'-g-' 

-[(1.65±ai6)s-^-g-2]ms, 

and 

/?« = (85.415±0.034)s-^-g-^ 

+ [(0.018±0.057)s-^-g-2]ms 

with V = (60-2) degrees of freedom. The mean /?„ 
values across all measurements and samples are 
/?a(Beckman) = (85.492 ±0.212) s"^ • g"^ and 
/?a(Packard) = (85.424±0.124) s"' • g'^ Several 
distinct observations may be noted. First, the LS 
response for the ^^Pb ce (as seen by the Rac-Ra 
differences) is greater for the Packard instrument, 
and is more dependent on changes in sample 
quenching (as seen by the magnitudes of the fitted 
slopes dRaJdnts). Second, the magnitudes of dRJ 
dms with attendant uncertainties indicate that /?« is 
independent of rris (and hence quenching). Third, 
all of the above four /?„ values, whether obtained 
from the regression intercepts or from the averag- 
ing, are statistically equivalent. 



2.5 1994 4170 LS Measurement Results 

The results of 1994 LS measurements, obtained 
with samples PI through P5, Ql through Q4, and 
Q5 through Q8 (see Table 1) were essentially a 
redundant and independent reconfirmation of the 
reanalyzed 1993 results. The LS samples had sys- 
tematically varied sample quenching using differ- 
ent sample compositions. As indicated in Table 1, 
some of the sample series used larger scintillator 
volumes to further minimize quenching. Analysis 
procedures used to obtain the counting rate con- 
centrations Rac and Ra were identical to those em- 
ployed with the 1993 data. The LS spectra 
integrating regions AEa (in terms of approximate 
beta energies) were: 

for the Beckman instrument 

^Ea^32 keV to 4400 keV (P samples) 

25 keV to 4400 keV (Q samples) 

and for the Packard instrument 

19 keV to 2000 keV (P samples) 

25 keV to 2000 keV (Q samples). 

Each sample was measured with each instrument 
from 5 to 18 times. All results were corrected for 
decay (as well as for the ^^^Po impurity) to the com- 
mon reference time of 1200 EST 15 March 1994. 

2.5,1 P and Q Series Samples Figure 13 sum- 
marizes the four sets of results obtained with the 
two LS systems. As previously: (1) Rm is invariably 
larger for the Packard instrument than that ob- 
tained with the Beckman (which had a longer pulse 
resolving time); (2) variations in Rac with sample 
quenching (reflecting differences in detection and 
resolution of the ^^^Pb ce) is also greater with the 
Packard; (3) values of Ra obtained with either in- 
strument are invariant of sample quenching; and 
(4) the Ra values for both instruments are statisti- 
cally equivalent and virtually identical. 

Table 4 contains a summary of the various means 
for Rae and Ra, along with their precision estima- 
tors, and comparisons of the fitted parameters for 
the unweighted linear regressions shown in Fig. 13. 
In these cases, the regressions were performed us- 
ing the mean counting rate concentrations for each 
sample. The uncertainties in the fitted parameters 
thus have v = (5 - 2) and v = (4 - 2) degrees of free- 
dom for the P and Q sample series, respectively. 

2.6 LS Spectral Shape Considerations 

Extensive examinations of the LS spectra were 
performed throughout various phases of this work. 
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Fig, 13. LS counting rate concentrations R^t, and R^ obtained with the two LS systems for the P and 
Q series samples in 1994. Closed squares (R^) and open squares (7?„) represent the mean values for 
samples Q5 through Q8 with the Packard; closed and open triangles represent J?„c and R^^ respec- 
tively, for samples PI through P5 with the Packard; closed and open triangles (R^ and R^) are for 
samples Ql through Q4 with the Beckman; and closed and open circles (Rm and R^) are for samples 
PI through P5 with the Beckman. Each plotted value corresponds to the mean of 5 to 18 replicate 
measurements on each sample. The error bars represent standard deviation uncertainty intervals on 
the means. The solid lines are unweighted linear fits to the data. Although the Rac values vary with 
the instrument used to perform the measurements (Packard or Beckman) and with sample composi- 
tions, all of the Ra values are statistically equivalent and invariant. 

Table 4. LS measurement results for the P and Q sample series obtained with the Beckman and Packard counting 
systems 







Beckman 






Packard 






Samples P1-P5 


Samples 


Q1-Q4 


Samples 


P1-P5 


Samples 


Q^Q8 


Parameter 


^«e' 


RJ" 


R.c' 


RJ" 


^«/ 


^«^ 


^«c« 


R..' 


«/ 


5 


5 


4 


4 


5 


5 


4 


4 


"m"' 


62 


62 


48 


48 


25 


25 


71 


71 


Mean R' 


87.787 


85.532 


88.011 


85.460 


89.278 


85.553 


89.847 


85.528 


(s-'-g-) 


















s' 


0.226 


0.175 


0.512 


0.310 


0.323 


0.264 


0.540 


0.184 


(s-'-g-') 


















Sm» 


0.033 


0.026 


0.084 


0.052 


0.072 


0.062 


0.071 


0.026 


(%) 


















Intercept iJ*" 


87.967 


85.553 


88.301 


85.374 


89.508 


85.536 


90.643 


85.518 


(s-.g-) 


















V 


0.059 


0.040 


0.254 


0.057 


0.047 


0.042 


0.117 


0.043 


(s-^.g-0 



















* Counting rate concentration obtained from integrating the full-energy LS spectra (see Sec. 2.4.1), with ^'***Po and decay 
corrections to the reference time 1200 EST 15 March 1994. 

b 209pQ ^ counting rate concentration (see Sec. 2.4.1), with ^****Po and decay corrections to the reference time 1200 EST 15 

March 1994. 

*^ Number of samples in the sample series. 

^ Total number of measurements on all samples in the series. 

' Mean counting rate concentration. 

^ Standard deviation in the counting rate concentration for the n determinations; v = (n^ - 1) degrees of freedom. 

« Relative standard deviation of the mean in percent for the n determinations; Jm = 100 s/m Vrtm- 

^ Fitted intercept value (at m^ = 0) for the counting rate concentration from unweighted linear regressions on the mean R 

for each of the n^ samples in the series (see Fig. 13). 

* Standard deviation for the fitted intercept value; v = (/Js - 2) degrees of freedom. 
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These examinations involved systematic studies of 
the detailed shapes and structure of the spectra as 
functions of the LS sample compositions, chemical 
quenching changes, and source configurations. 
They were performed largely in order to under- 
stand the observed LS measurement results (par- 
ticularly in terms of interpretations and 
compatibility with the ^Po decay scheme), and to 
insure the validity of the analysis procedures used 
to obtain the previously given Ra values. 

2.6,1 Quench-Dependent Gross Effects Typi- 
cal gross effects of chemical quenching on the ^^Po 
spectrum are illustrated in Fig. 14 for the Q5 
through Q8 samples. The four spectra were ob- 
tained with the Packard instrument. The quench 
parameters varied from tSIE =528 for least- 
quenched sample Q5 to tSIE =414 for Q8 and 
were determined with a precision having a relative 
standard deviation of the mean of 0.16 % to 
0.17 %. Two dominant features can be noted. With 
increasing quenching in the samples (i.e., decreas- 
ing tSIE values), the broad a peaks at around an 
equivalent beta energy of 350 keV to 550 keV shift 
to lower energies, and the integrated areas of the 
low-energy ce peaks progressively decrease in mag- 
nitude. These latter peak areas are, of course, 
merely the Re values as given, for example, by the 
changing (Rac^Ra) differences in Fig. 13. The a 
peak areas (given by/?„ in Fig. 13), although shifted 
to different energies, are invariant of the quench- 
ing changes. 

These a peak shapes are considered in greater 
detail in Fig. 15. One may observe that the spec- 
trum for the least quenched Q5 sample almost has 
the appearance of an unresolved doublet. This fea- 
ture was also seen (at times in this laboratory, as 
well as by other researchers) in other LS spectra 
from monoenergetic a-emitting radionuclides, e.g., 
^"**Po [3], and is frequently attributed to "optical 
effects" arising from refractions of the light pho- 
tons in the LS vial walls, or from light reflections in 
the LS counting chamber or off the vial walls. Simi- 
lar "doublets" were also observed however, in the 
earlier 1993 spectra obtained with polyethylene 
vials, and was seen in spectra from both the Beck- 
man and Packard LS counting systems. These were 
surprising findings. One would think that the glass 
and polyethylene vials have considerably different 
refraction and surface reflection properties; and 
that there undoubtedly are some differences in the 
"optics" for the two counting systems. This "dou- 
blet" effect is rarely apparent in more highly 
quenched samples when the alpha peak shifts to 



lower energies as shown in Fig. 15.^ With the shifts 
to lower energies, the resolution (or peak widths) 
also decrease on an absolute basis. The peak 
widths, given by full-width-at-half-maximum 
(FWHM) values, for the four spectra of Fig. 15 
range from about 120 keV to 170 keV. On a rela- 
tive basis compared to the peak centroids, how- 
ever, all four have FWHM values of approximately 
30 % (within 1 % or 2 %). 

2.6.2 ^^Anomalous Bumps" Another surpris- 
ing feature of the LS spectra is the presence of 
"anomalous bumps" distinctly separated from the 
a peaks and at energies just below these peaks. 
These bumps were actually identified prior to un- 
derstanding the detection of the delayed 2.3 keV ce 
when the LS results still seemed inexplicably dis- 
crepant. It was thought that perhaps the differ- 
ences in the initially observed R^ counting rates 
between the two LS counters might be due to these 
mysterious bumps that could not be attributed to 
any known transitions in ^^Po decay. As a result, 
the spectral structures and magnitudes of these 
bumps were extensively examined and systemati- 
cally studied. They were observed in all of the LS 
spectra, in all of the samples ranging from the 1993 
A ,B, M, and N series of samples to the 1994 P and 
Q sample series (which therefore encompassed use 
of both polyethylene and glass LS vials), and in 
both counting systems. 

Figure 16 contains a typical example of these 
"anomalous bumps" as seen in spectra of sample 
P3 obtained with both spectrometers. The bumps, 
like alpha peaks, shift to lower energies with in- 
creased quenching as seen in the series of spectra 
in Fig. 17. In this case, the bump centroids decrease 



^ Other explanations can not be definitively excluded. In a pri- 
vate communication, Dodson [17] reported that he has only ob- 
served this type of doublet effect with ^'^^Am, and not with other 
a emitters like ^^Ra or ^^Pu. He is lead to believe it may be a 
sample preparation effect since he never observed it with sam- 
ples prepared by extractions (i.e., with extractive scintillators). 
In a set of quenched ^*^Am "standards", a doublet was "not 
initially present, but later the doublet appeared at the same 
time that the emulsifier precipitated, similar to a two-phase con- 
dition" [17]. This explanation, however, is not compelling for 
the present situation. Firstly, the observed doublets in the ^"^Po 
spectra appear to be composed of two symmetric components of 
nearly equal magnitude. It is hard to imagine that the polonium 
was in different chemical states or phases to nearly exactly equal 
extents. Secondly, the observed doublets were most pronounced 
in the least aged and least quenched samples, and not the con- 
verse. Although the "doublets" were observed in spectra ob- 
tained with both spectrometers, their appearance was invariably 
more pronounced (or defined) in those obtained with the Pack- 
ard instrument. 
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Fig. 15. Details of the broad alpha peaks shown in the full spectra of Fig. 14. The peak widths (FWHM) 
on a relative basis and peak areas are approximately equal in all four samples Q5 through Q8. 



from about 210 keV in the least quenched sample 
Q5 to about 150 keV in sample Q8. The data of 
Fig. 17 (unlike all other displayed LS spectra) were 
smoothed using an averaging Savitsky-Golay differ- 
entiation technique to eliminate the wide fluctua- 
tions in the data (arising from large variations in 
the small counting rates per channel and their in- 
herently large statistical "counting errors") and 



thereby to more clearly show the bump structures. 
The shifts of the bumps to lower energies were not 
found to be systematically correlated to the a peak 
shifts. About the only conclusive statement that 
can be made is that the energy differences in the 
bump shifts were less than that for the shifts in the 
a peak centroids for given quench changes. If a 
correlation exists, it is not a simple relationship. 
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Fig. 16. Comparison of the "anomalous bumps" in the spectra 
for sample P3 obtained with the two LS counting systems. 




Approximate Beta Energy (keV) 

Fig. 17. Comparison of the "anomalous bumps" in the LS spec- 
tra of four samples having different quench conditions. 

There were also suggestions that the bumps be- 
came more dispersed and broadened (on both an 
absolute and relative basis) with increasing sample 
quenching, and that perhaps even secondary 
bumps appeared in the valleys between the princi- 
pal bump and the a peak (see, for example, the 
spectra for Q7 and Q8 in Fig. 17). The relative 
magnitudes of the bumps in comparison with the a 
peak (in the same spectrum) were highly variable, 
ranging from less than 0.32 % to over 2.2 %. Again, 
there was no obvious correlation or simple rela- 
tionship between the magnitude of the bumps and 
any sample quenching or sample counting condi- 
tions. 



For reasons outlined immediately below, it is be- 
lieved that these anomalous bumps are most likely 
due to optical effects from a interactions in the 
cocktails, and not due to LS responses from some 
other non-a transitions. Why, how, when, and 
where they appear in the spectra, however, still re- 
mains unknown. Sample preparation effects due to 
polonium in a two-phase condition can not be 
ruled out as the causal factor. This would be sur- 
prising however considering the presence of the 
bumps over such a wide range of sample composi- 
tions and conditions. It should also be emphasized 
that if these bumps were not included as part of 
the 4iTa LS response, then the /?« values given pre- 
viously would not have formed a quench-indepen- 
dent and consistent data set. 

2.6.3 Alpha-Beta Pulse Discrimination Sev- 
eral LS spectra obtained using alpha-beta (a-p) 
pulse discrimination were also measured to insure 
that the previously ascribed spectral interpretations 
were valid, particularly in regard to the assign- 
ments of the "anomalous bumps" to a interactions 
and of the very low ~3 keV peak to conversion 
electrons. 

The LS counting sources used to obtain these 
spectra consisted of about a 1 g aliquant of the 
^^Po standard solution (obtained from a randomly 
selected sealed ampoule, #146) in about 14 g of a 
prepared scintillation fluid. The scintillating cock- 
tail was "Ultima Gold AB" (Packard Instrument 
Co.) which was selected because it provides slow 
pulse decay characteristics which are necessary for 
effective a-P pulse discriminations. The prepared 
scintillator is also reported to be very efficacious 
for accepting high mineral acid loadings (up to 2 
mL of 2 mol • L'^ HNO3 in 10 mL of scintillator). It 
may be noted that the previously used "Ready 
Safe" scintillation cocktail did not provide ade- 
quate a-P pulse discrimination. 

The spectra were obtained with a slightly differ- 
ent model of the Packard Tri-Carb 2500TR LS 
spectrometer that was used for all of the previously 
described Packard instrument measurements. 
Pulse discrimination between a and p pulses relies 
upon their differing pulse decay times. Initial spec- 
tra were taken with a nearly arbitrarily chosen mid- 
range pulse discrimination setting with the 
expectation that there would be some accidental 
overflows of a pulses in the p spectrum, and vice 
versa. 

Figure 18 gives a composite of the respective a 
spectrum (broken line) and p spectrum (solid line) 
obtained with the initially chosen discriminator set- 
ting. As shown, the low energy peak centered at 
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Fig. 18. LS Spectrum of ^"^Po obtained using alpha and beta 
pulse discrimination, 

around 3 keV is clearly due to 3 (i.e., electron) 
pulses; whereas a selected (and 3 discriminated) 
pulses clearly comprise the ascribed broad a peak. 
Expanded scale details of these spectra in the en- 
ergy interval of 50 keV to 600 keV are given in the 
upper trace of Fig. 19. Of the intense a peak, only 
a very small fraction of the response is from 3 se- 
lected pulses. These are attributed to "accidentals" 
whose relative magnitude was dependent on the 
specific discrimination setting used. At this same 
setting, less than 20 % of the total area of the 
"anomalous bump" was due to selected 3 pulses. 
Spectra obtained with a 25 % change in the initial 
a-3 pulse decay discrimination setting resulted in 
those in the lower trace of Fig. 19. In this latter 
case, the accidentals under the a peak dramatically 



decreased and virtually the entire anomalous bump 
is comprised of selected a pulses. The findings re- 
inforced the belief that the anomalous bumps were 
indeed the result of some unknown optical effects 
in detecting some small fraction of the a interac- 
tions. 

2.6.4 Optical Effects in DifTering Counting 
Source Configurations Lastly, in an attempt to 
discern (if not understand) the possible optical ef- 
fects, several spectra were obtained using sources 
of very different counting configuration. It was 
thought that the "optics" of these configurations 
would of necessity be substantially different. The 
configurations consisted of a normally prepared 
sample source containing mtot=*ll g of scintillator 
and sample in a 22 mL glass LS vial (nominal 2.5 
cm diam.); a source containing mxoi—23 g of scintil- 
lator and sample in a similar 22 mL glass vial such 
that the vial was completely filled with cocktail; 
and another completely filled source in a smaller 
1.4 cm diameter glass vial that was in turn concen- 
trically contained in a CeHu (hexane) filled 22 mL 
vial. The sources are designated N (for "normal"), 
F (for "full"), and D (for "double"), respectively. 
The sources were prepared with closely matched 
chemical quenching conditions, and the spectra ob- 
tained with them were normalized to the ^**^Po solu- 
tion mass contained in each. In this way, if the 
sources had been in identical configurations and 
with identical optics, then one could expect identi- 
cal spectra for them. 
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Fig. 19. Details of the LS spectrum of Fig. 18 showing the relative alpha and beta re- 
sponses for the "anomalous bump" and a peak at two pulse discrimination settings. 
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Comparison of the N and D sources was made to 
see if the **anomalous bumps" were at all associ- 
ated with the presence of an air space in the N 
source configuration. Conceptually, it was thought 
that perhaps a lower energy a peak might arise be- 
cause of scintillations from any cocktail (above the 
air gap) that adheres to the vial cap surfaces. The 
filled configuration of source D would obviate this 
condition, although there would of necessity be 
some loss of light from scintillations occurring in 
the cocktail volume above the vial cap. Compari- 
sons were made with the D source configuration 
since it is known that the resolution of a spectra 
can be improved by decreasing the sample source 
cross-sectional area (e.g., by using smaller diameter 
sample tubes) [22]. The smaller diameter vial of 
source D was placed inside the larger 22 mL vial to 
keep the total source volume constant. The glass 
wall of the inner vial was virtually transparent since 
it was surrounded by hexane, which has a nearly 
identical refraction index as the scintillation cock- 
tail. 

The results of these comparisons are shown in 
the spectra of Figs. 20, 21, and 22. The first (Fig. 
20) provides the full spectra including the low en- 
ergy ce peaks. The second (Fig. 21) gives details of 
the a peaks. The broad a peak tailing into the 
"anomalous bump" region is shown lastly (Fig. 22). 
The integral intensity of the ce peaks progressively 
decreases in going from the N to F to D configura- 
tions. This may partly be due to light losses above 
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Fig. 20. Comparative ^^^''Po LS spectra for three very different 
source configurations N, F, and D. Refer to text for details. The 
upper trace compares the spectra for N and F; the lower com- 
pares N and D. 
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Fig. 21. Details of the broad alpha peaks for the three source 
configurations shown in the spectra of Fig. 20. 
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Fig. 22. Details of the "anomalous bumps" for the three source 
configurations shown in the spectra of Fig. 20. 

the vial caps in the fully filled F and D configura- 
tions. The N configuration shows the apparent 
"doublet" nature of the broad a peak and the by- 
now consistent "anomalous bump.'* In source con- 
figuration F, the bump is still present but is of less 
intensity and is shifted to a lower energy. At the 
same time, the a peak from source F has substan- 
tially broadened in comparison to source N, and 
has even more of the appearance of an unresolved 
doublet with a very long low energy tail. The most 
surprising results were that for source configura- 
tion D. The a resolution with a peak centered at 
350 keV was as expected considerably improved 
with a FWHM of about 70 keV which is less than 
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half that obtained with configuration N. Yet, this a 
peak is also accompanied by another large and ad- 
joint secondary a peak centered at 250 keV. This 
secondary peak encompasses a large fraction of the 
total a response (about 1/4) and also contained 
very long low energy tailing. The spectra of config- 
uration D did not exhibit any of the usually-found 
"anomalous bumps," but the tailing on the sec- 
ondary a peak had a discrete inflection point. All 
of the above described spectral structures (in every 
source configuration) appear, in fact, to be nothing 
but a complex variety of shape changes (of a con- 
tinuous nature) in the a peak. 

These spectral shape findings can not be ex- 
plained (at least not by the present authors) in 
terms of any overall qualitative theory or model for 
the scintillation optics. The findings, however, 
tended to support previous contentions that the 
"anomalous bumps" were associated with true a 
interactions in the scintillator, and that they arise 
from some type of light-collection optical effects 
that in turn could be influenced by changes in 
source configuration. 

2.7 Confirmatory Measurements by lira 
Proportional Counting 

Independent confirmatory measurements of the 
^Po a emission rate concentration were per- 
formed by lira proportional counting of solid point 
sources LXl and LX2. The detector consisted of a 
12 cm diameter hemispherical gas-flow propor- 
tional counter. The counting gas was a "P-10" 
argon-methane (90:10) mixture. A detailed de- 
scription of the detector and its operation has been 
previously described by Hutchinson [23]. 

The results of duplicate measurements of each 
source are given in Fig. 23. The illustrated a emis- 
sion rate concentrations as a function of source 
mass were obtained from the observed net count- 
ing rates (after background subtraction) which in- 
cluded an extrapolation to zero pulse height [23]. 
The corrected net rates were converted to concen- 
trations by applying the known solution masses that 
were used to prepare the sources. The detection 
efficiency was assumed to be exactly 0.5 counts per 
alpha for a perfect 2t: steradian solid angle, 
exclusive of source self-absorption losses and 
source-mount back-scattering effects. The relative 
standard uncertainty in the 2iv efficiency is typically 
taken to be ±2.5 % [23]. The sources had been 
measured in mid-June 1993, but the values re- 
ported here were decay corrected to the 15 March 
1994 reference time for direct comparison with the 
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Fig. 23. Confirmatory measurements of the ^^'^Po a-emission 
concentration by lira gas-flow proportional counting. An ex- 
trapolation to zero source mass (without backscattering correc- 
tions) is given by a linear fit to the data (solid line). The broken 
lines are the uncertainty limits for the regression. 

LS measurement results. These values also contain 
a small correction for the known ^*^Po impurity (in- 
cluding its decay) so that they correspond to only 
the ^^Po a emission rate concentration. The error 
bars shown on the data of Fig. 23 represent one 
standard deviation uncertainty intervals for the as- 
sumed Poisson-distributed total statistical "count- 
ing error." These uncertainties ranged from 
±0.3 % to ±1.1 % for the four determinations. 

Figure 23 also contains an extrapolation of the a 
emission rate concentration to zero source mass to 
account for source self absorption. The solid line 
corresponds to an unweighted ;^^-minimized linear 
regression fitted to the data. The uncertainty 
bounds for the fit are shown as dotted lines. The 
fitted intercept (at zero source mass) for the a 
emission rate concentration is (88.41 ±0.56 s"^ • 
g"V Using a weighted linear regression (with 
weighting by the reciprocal of the Poisson "count- 
ing error" variance), the fitted intercept is 
(88.52 ±0.52) s"^ • g"^ which is indistinct from that 
for the unweighted fit. For "weightless" sources 
mounted on platinum, the backscattered fraction in 
2iT has been shown to be 0.0175 ±0.0025 [24,25]. 
Invocation of this 3.5 % back-scattering correction 
results in a corrected value of 85.5 s"^ • g~^ for the 
^'^^Po a emission rate concentration for the solution 
standard. These applied corrections for self ab- 
sorption (by extrapolation to zero mass) and back 
scattering (which assumed a condition of source 
weightlessness) are of course somewhat simplistic. 



^ In correct usage, alpha emission rate concentrations have units 
of s"^ • g~*, rather than the traditional units of "alphas per 
second per gram." 
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In actuality, the source mass extrapolation of Fig. 
23 included some part of backscattering effects 
since the magnitude of any source mount scattering 
will of necessity by dependent on the source mass 
itself. The two effects are complexly related. Nev- 
ertheless, the relative standard uncertainty in the a 
emission rate concentration may be considered to 
be roughly about ±3 %, considering the uncertain- 
ties in the fitted intercept and in the 2'n-a detection 
efficiency (including the source absorption and 
scattering effects). Clearly, the 2'rra measurements, 
although less accurate, admirably served as a suit- 
able confirmatory determination that was com- 
pletely independent of the LS measurement 
results. 
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Fig. 24. Comparative a spectra for sources LXl and LX2 ob- 
tained with a Si surface-barrier junction detector. The spectra 
were normalized to the same peak heights to show the greater 
low-energy tailing in the more massive LXl source. 



2.8 Alpha-Emission Impurity Analyses 

Analyses for a-emitting radionuclidic impurities 
were also made with the LXl and LX2 point 
sources. The analyses were performed by a spec- 
trometry using a Si surface-barrier junction detec- 
tor which was suitably coupled to a 1024-channel 
pulse-height analyzer. Conversion gain settings 
were typically 5 keV to 10 keV a energy per chan- 
nel. The sources were measured through a very 
thin film absorber to prevent surface contamina- 
tion of the detector face. The film consisted of 
VYNS, a polyvinylchloride— polyvinylacetate co- 
polymer, of nominal surface density « 10 |xg • cm"^ 
that was stretched on an aluminum annular sup- 
port ring. Multiple spectra with counting time in- 
tervals ranging from 2x10"* s to 6x10^ s were 
obtained over the course of 26 d in the time period 
of May to June 1993. 

Typical spectra for the two sources, showing the 
partial region for the principal ^Po 4.883 MeV a 
transition, are given in Fig. 24. The effective energy 
resolution for the detector/source combination, in 
terms of the peak full width at half maximum 
(FWHM), is roughly 50 keV to 60 keV. This can be 
contrasted to an optimum 10 keV to 30 keV which 
would be expected from a "weightless" source 
without a backing mount. Low energy tailing, re- 
sulting from degradation of the a full energy as 
well as from some possible small influence of the 
unresolved ^Po 4.622 MeV a, is evident in the 
spectra from both sources. The weak 4.622 MeV a 
[/a = (0.0092 ±0.0005) alphas per decay] is com- 
pletely masked by the more dominant 4.883 MeV 
peak. The low-energy peak broadening and ob- 
served FWHM energy resolution can not be at- 
tributed solely to ^Po source material self 
absorption since the solution deposits used to pre- 



pare the sources were carrier free and therefore 
should have been nearly "weightless" after drying. 
More probable contributing factors were: inhomo- 
geneous source depositing; deposits of unknown 
chemical impurities (most likely Pb and Bi); 
backscattering from the thick Pt source mounts; 
surface irregularities in the back mounts; and ab- 
sorption in the VYNS barrier. 

The only observed a-emitting impurity, based on 
a total cumulative counting time exceeding 
1.6x10^ s (-19 d), was ^'^^'Po whose principal a 
(/„ = 0.999982 ± 0.000002) has an energy of 
£a = (5.1152±0.0015) MeV [11]. There was no evi- 
dence of ''"Po presence [£„ = (5.30438 ±0.00007) 
MeV; /a = 1.00] [26,27] which is unsurprising be- 
cause of its substantially shorter half-life (r = 138.4 
d) and age of the ^**^Po stock material. Detection of 
these two most likely polonium isotope impurities 
is greatly facilitated in that their principal a transi- 
tions have energies greater than that for ^^o. Fig- 
ures 25 and 26 illustrate the relative ^Po and ^^Po 
a intensities. The former, in full scale, clearly 
shows the large relative magnitude difference in 
the two peaks, although that for ^***Po is barely per- 
ceptible; the latter, in an expanded scale, demon- 
strates the distinct resolution and definition of the 
5.115 MeV a from the ^Po impurity. 

Three independent determinations of the ^"**Po 
to ^**^Po impurity ratio, in terms of their respective 
a-particle emission rate concentrations, were per- 
formed. The results, which were decay corrected to 
the common reference time of 1200 EST 15 March 
1994, are summarized in Table 5. The propagated 
relative uncertainties for the "counting errors" and 
for location of peak backgrounds ranged from 
±1.6 % to ±3.9 % and from ±1.6 % to ±2.1 %, 
respectively. The ^Po/^**^Po impurity ratio mean 
for the three determinations was 0.00124 with a 
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Fig. 25. Detail of the ^^Po a spectrum showing the relative magnitude of the ^'*'Po impurity. The spectrum was 
obtained with source LX2 and the Si surface-barrier junction detector. 

MeV, Xj^Fo^ 6 x 10"^; and in the region E^^ 3.5 
MeV, X«^'Po^2xlO-^ 



2.9 Photon-Emission Impurity Analyses 

Analyses for photon-emitting radionuclidic im- 
purities were performed with samples Gl through 
G3 using two different spectrometry systems. The 
systems comprised (a) an intrinsic p-type Ge detec- 
tor (having a Nal(Tl) equivalent efficiency of 18 %) 
coupled to a 16384-channel pulse height analyzer 
(PHA) that spanned a 40 keV to 3300 keV energy 
region; and (b) an intrinsic n-type Ge detector 
(12 % Nal(Tl) equivalent efficiency) which used a 
8192-channel PHA spanning 8 keV to 915 keV. 
Both systems employed a 60 Hz pulser for rate nor- 
malizations. Detection efficiency curves for the sys- 
tems were determined with calibrated (NIST 
SRM) 2^* Am, "^In, ^'Co, ^^Cr, «^Sr, *^'Cs, and ^Co 
sources. 

Photons associated with both the EC and a 
branch decays of ^^Po include Bi K x rays arising 
from K capture in the EC branch and from internal 
conversion fluorescence of the 896.6 keV transition 
in ^Bi; the 896.6 keV 7 ray in ^^Bi; Pb K x rays 
from internal conversion of the 260.5 key and 
262.8 keV 7 rays in ^^Pb; and the 260.5 keV and 
262.8 keV 7 rays in ^^Pb (see Fig. 2). All of these 
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Fig. 26. Expanded scale detail of 5.115 MeV a from the ^"*'Po 
impurity in the ^^Po a spectrum of Fig. 25. 

standard deviation of 0.00010 (or ±8.1 % on a rel- 
ative basis). 

Further examinations of the spectral back- 
grounds in certain energy regions could establish 
limits for the presence of any other a emissions. 
Impurity ratios, in terms of emission rate concen- 
trations Xj^^Fo, for a particles of energy £„ (aris- 
ing from an a-emitting impurity Xa in the ^Po 
standard solution) have the following lower limits: 
in the region £:a>5.18 MeV (which includes ^^^Po), 
Xc/**^Po< 2x 10-^ in the region 3.5 M6V^Ea^ 4.2 
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Table 5. 


Determinations of the ^'^PoF^Fo 


impurity ratio by a 


spectrometry 






Source 
identity 


Counting 
time 


209po 

4.88 MeV a" 
(s-i.g-0 


208po 

5.11 MeV a^ 
(s-^.g-0 


% 


% 


208po;209poa 

Impurity ratio 


LX2 
LX2 

LXl 


3.5 X 10^ 
6-0 xlO^ 
6.0 xlO^ 


22.40 
22.07 
21.46 


0.0256 
0-0297 
0.0265 


3.9 
2.8 
1.6 


2.1 
1.6 
2.0 


0.001144 
0-001344 
0.001236 



■ As of the reference time 1200 EST 15 March 1994. The counting rate concentrations for the ^**^Po and 
^^^''Po peaks are in traditional units of "counts per second per gram." 

^ Relative standard deviation in percent for the assumed Poisson-distributed tot al statistical "counting 
error" which is given by Sp - 100 V(^ai + Nbi)/{Nai - Nbif + (Naz + A^b2)/(^«2 - N^zf where N^i and N^z are 
the total gross counts in the two a peaks and M>i and M>2 are the estimated total counts in the two 
underlying and extrapolated peak backgrounds. 

*= Relative standard uncertainty in percent for locating and determining the two underlying a peak back- 
grounds (propagated from estimated uncertainties on each peak). 



photons, including partial resolutions of the Pb and 
Bi Kaia2 and Kpip2 x rays, were detected. Determi- 
nations of their absolute intensities, based on the 
present ^^Po calibration, will be reported sepa- 
rately [28] in a re-evaluation of the ^^Fo decay 
scheme. 

No photons other than those associated with the 
decay of ^**^Po were detected. Emission rate con- 
centrations^ for any other photons arising from 
photon-emitting impurities in the ^Po standard 
solution have the following lower limits: < 2 x lO"'* 
s~^ • g"^ in the energy region 15 keV to 68 keV; 
< 1.5x10"'* s"^ » g~^ in the region 81 keV to 256 
keV; <6x 10"^ s~^ • g"^ in the region 266 keV to 
892 keV to; and <4xlO~^ s"^ • g"Mn the region 
900 keV to 3300 keV, provided that the energy of 
any such unobserved photon is 4 keV or more dif- 
ferent from the emissions accompanying the ^°^Po 
decay. These limits were based on numerous inde- 
pendent measurements of the sources with the two 
spectrometer systems over a total counting time in- 
terval of 7.2 x 10^ s (-83 d). 



2.10 Calibration Summary and Uncertainty 
Analysis 

2.10.1 Mean LS Counting Rate Concentration 

The LS determinations of Ra obtained with both 
counting systems and with all of the 1993 and 1994 
samples are summarized in Table 6. The means 
given in Table 6 were obtained from unweighted 
averages across all measurements in the given sam- 
ple series. Weighted means were generally indistin- 
guishable since reciprocal variance weighting 
factors were all of the same nominal magnitude. 



^ These emission rate concentrations have traditional units of 
''photons per second per gram." 



As shown previously in Figs. 8 through 13 and 
Tables 2 through 4, the mean value of Ra for any 
sample within a given sample series and obtained 
with a given instrument was statistically equivalent 
(with the possible exception of sample B6— see 
Figs. 8 and 9) with that for any other sample. Simi- 
larly, as shown in Table 6, when these sample 
means are averaged over a given sample series, the 
resulting means for any two sample series were also 
equivalent. Further, the means obtained with ei- 
ther the Beckman instrument or the Packard were 
invariant. 

The ratio of mean Ra for the A and B sample 
series to the mean Ra for the M through Q sample 
series may be directly compared to the gravimetric 
dilution factor. These are: 

D (Beckmanl- (12560±27) s"^ ■ g"^ 
D (Beckman) -^g^^g^^p 244) s"^ • g"^ 

= 146.95 ±0.52 

D (Packard)- (12553 ± 26) s"^ ■ g'^ 
V (FacKard)-^g5^gg^p 209) s"^ • g"^ 

= 146.88 ±0.47 

compared to D (gravimetric) = 147.24 ±0.07 which 
is approximately 0.2 % larger. Although they are in 
statistical agreement, the absolute differences are 
slightly larger than that usually obtained by our 
laboratory for verifications of gravimetric dilution 
factors. Typically, LS verification measurements of 
dilution factors are <0.1 %. As noted previously, 
the differences may in part be due to chemical 
quenching differences arising from greater quanti- 
ties of HNO3 in the A and B sample series. Never- 
theless, these observed differences are irrelevant 



27 



Volume 100, Number 1, January-February 1995 

Journal of Research of the National Institute of Standards and Technology 



Table 6. Summary of all LS 


measurement results for 


mean R, 


a obtained with Beckman and Packard i 


counting systems 














Results with Beckman 




Results with Packard 




Sample 


209po 




Number of 


n" 


Ra mean*^ 


^m 


n* 


Ra mean** 


^m 


identities 


Solution 


samples 




(s-^.g-0 


%' 




(s-^-g-0 


%<= 


A1-A2 


A 




2 


59 


12566 


0.016 


34 


12557 


0.019 


Bl-Bbx 


B 




9 


138 


12557 


0.020 


115 


12552 


0.021 


all A & B 


A/B 




11 


197 


12560 


0.015 


149 


12553 


0.017 


MA1-MA3 


M (amp. 


39) 


3 


30 


85.494 


0.057 


30 


85,418 


0,043 


MB1-MB3 


M (amp. 


93) 


3 


30 


85.448 


0.056 


30 


85.435 


0.051 


MC1-MC3 


M (amp. 129) 


3 


30 


85.333 


0.049 


30 


85.438 


0.043 


MD1-MD3 


M (amp. 205) 


3 


30 


85.464 


0.047 


30 


85.439 


0.045 


N1-N6 


M (amp. 


47) 


6 


30 


85.492 


0.045 


60 


85.424 


0.019 


all M & N 


M 




18 


150 


85.446 


0.023 


180 


85.430 


0.018 


P1-P5 


M (amp. 


40) 


5 


62 


85.532 


0.026 


25 


85.553 


0.062 


Q1-Q4 


M (amp. 


29) 


4 


48 


85.460 


0.052 








Q5-Q8 


M (amp. 


29) 


4 








71 


85.528 


0.026 


all P & Q 


M 




9 


110 


85.501 


0.027 


96 


85.535 


0.025 


AllM-Q 


M 




27 


260 


85.469 


0.018 


276 


85.466 


0.015 



* Total number of measurements on all samples in the series. 

^ Mean of the ^^Po a counting rate concentration, with ^^'Po impurity and decay corrections to the reference time 1200 EST 15 March 

1994 (see Sec. 2.4.1). 

^ Relative standard deviation of the mean in percent for the n determinations. 



insomuch as the ^**^Po cahbration can be based di- 
rectly on the measurements of samples (sealed am- 
poules) of the master solution. 

The overall results of Table 6 may be averaged in 
several ways. Consideration of the possible mean 
values for Ra derived in these different ways is in- 
formative. With the Beckman instrument, the 
mean Ra for all 27 samples of the M through Q 
sample series, involving a total of 260 independent 
LS determinations of /?« is (85.469 ±0.244) s"^ -g'^ 
The above standard deviation has a corresponding 
relative standard deviation of the mean of 
Sm= ±0.018 %. For the Packard instrument, based 
on 276 measurements with 27 samples, the mean 
Ra is (85.466 ±0.209) s"^ • g"^ (5^ = 0.015 %). As 
noted, these aggregated means can be compared to 
that obtained from replicate measurements on any 
one sample (Tables 2 and 3; and Figs. 8 through 
13), as well as to the means obtained from all sam- 
ples within the various sample series (Tables 2 
through 4, and 6). They are all equivalent. One 
could alternatively average across the means ob- 
tained from the results from each of the seven am- 
poules to obtain grand means. These are (85.460 ± 
0.063) $-' • g-' (5m = 0.024 %) and (85.460 ±0.055) 
s~' • g"* (5m = 0.021 %) for the Beckman and Pack- 
ard results, respectively. The precision estimators 
on these grand means reflect the between-ampoule 
variability. Lastly, the overall mean obtained from 



all 536 determinations (from both LS instruments) 
on all 27 samples from the 7 different ampoules is 
/?„ = (85.468 ±0.227) s"^ • g"^ (^.=0.012 %). 

It is significant to realize that the 536 measure- 
ment results embodied in this last mean were ob- 
tained under a tremendous range of conditions. 
They were obtained: (i) with two different counting 
systems; (ii) using different AEa integrating re- 
gions; (iii) with different counting time intervals; 
(iv) from measurements separated by more than a 
year in time; (v) with background subtraction vari- 
abilities over this time (vi) with decay corrections 
for ^**^Po as well as for the ^^Po impurity over this 
same time; (vii) using counting sources in both 
glass and polyethylene LS vials; (viii) with sources 
prepared from 7 randomly selected ampoules of 
the ^**^Po standard solution; and (ix) with 27 sam- 
ples having wide variations in sample composition 
(including aliquant size of the ^Po solution and 
hence emission rates) and chemical quenching con- 
ditions. 

The results of these determinations of Ra were 
also subjected to extensive statistical tests on sub- 
sets of the data across a variety of variables: sample 
identity including sample preparation order, sam- 
ple composition, sample masses, sample replicates 
for a given ampoule of the ^^Po solution, the tim- 
ing and sequence of the LS measurements, and the 
LS measurement system used. These included 
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;^^-tests of the homogeneity in subsets of the ob- 
served sample variances (across the variables), F- 
tests of the homogeneity in the various subset 
sample means, r -tests of differences between the 
various subset means, and tests of possible correla- 
tions and biases using analysis of variance 
(ANOVA) techniques with sequential two-way 
classifications of the variables. None of the tests 
indicated any statistically significant differences in 
any of the tested subset sample means or sample 
variances. The agreement in the results between 
the various sample series and between the two LS 
counting systems, as provided in Table 6, are repre- 
sentative. 

The 536 values comprised in the mean Ra = 
(85.468 ±0.227) s'^ • g"^ are normally distributed. 
This is demonstrated in Fig. 27 where a frequency 
distribution of the data is superimposed on a Gaus- 
sian (normal) probability distribution function 
(solid line) having parameters of ^i = 85.468 and 
{r = 0.227.^ The probability function was normal- 
ized to the n =536 sample size so that both distri- 
butions have equal areas. 



,1* 
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Fig. 27. Frequency distribution for the 536 determinations of 
the LS counting rate concentration R^. The solid line is for a 
normal distribution having parameters of /i =85.462 s"' • g"^ 
and o- = 0.227s-' • g-^ 



** Several standard statistical tests were performed to test the 
data set for randomness and for its distributional form. These 
included both graphical evaluations and calculations of various 
distributional measures. The 536-value data set appeared very 
normally distributed on a normal probability plot, particularly if 
the two largest values (see Fig. 27) in the set were excluded. 
Tests for normality using a probability plot correlation coeffi- 
cient test gave coefficients of r = 0.998 for the slightly smaller 
n =534 data set and r =0.996 for the entire {n =536) data set. 



2.10.2 a Emission Rate Concentration The 

^Po a-particle emission rate concentration was ob- 
tained from the mean Ra on application of the two 
aforementioned correction factors, viz., Aa = 

The first, /:ze, was found to be highly variable 
from spectrum to spectrum because of the low 
counting statistics and imprecision in locating a 
baseline under the ce peak in extrapolating to zero 
energy. As a result of this variability, it was decided 
to apply an average correction to all of the spectra. 
Values of kzE ranged from lows of /:ze==' 1.0002 to 
highs of /:2E=^ 1.0008, with a majority centering 
around a value of /:ze=^1.0006 (i.e., the integral 
area for the extrapolated baseline under the ce 
peak was 0.06 % of Ra). No obvious correlation 
between the observed magnitudes of kzE and 
quenching could be found. The estimated standard 
uncertainty on the central value of 1.0006 was 
taken to be ±0.0003 (i.e., ±0.03 % on a relative 
basis). There is an additional uncertainty, esti- 
mated to be ±0.05 %, for extrapolation to zero 
pulse height. The relative combined standard un- 
certainty in kzE is thus taken to be ±0.06 %, such 
that feE= 1.0006 ±0.0006. 

The second correction factor, k^c = 0.9988 ± 
0.0006, is for the contribution to the LS response 
from radiations that arise from the ^Po EC decay. 
The above value was obtained from approximations 
given in Appendix A. 

The mean a emission rate concentration is there- 
fore y4„=85.42 s"^ • g-^ (as of the 1200 EST 15 
March 1994 reference time). Assuming an a-decay 
branching ratio of /„ = 0.9952 ±0.0004 from Martin 
[11], the corresponding ^*^^Po activity concentration 
is 85.83 Bq-g-\ 

2.10.3 Uncertainty Assessment The relative 
combined standard uncertainty in Aa is Wc = 
±0.17 %. The relative expanded uncertainty, using 
a coverage factor ofk=2,\sU-± 0.34 % which is 
assumed to provide an uncertainty interval of 90 % 
to 95 % confidence. The component uncertainties 
used to obtain these estimates are outlined in 
Table 7. These component uncertainties are given 
in terms of their contribution to the relative stan- 
dard uncertainty in ^a, i.e., UA„(j) = [(dAjdzj)u,fiyAa 
for each component Zj having a relative uncertainty 
estimate of Wr,. Of necessity, some of the uncer- 
tainty components are difficult to quantify on a 
purely objective basis and are therefore based on 
subjective experiential estimates. Analyses, justifi- 
cations, and underlying assumptions used to quan- 
tify all of the components follow. 
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Table 7. Summary of the uncertainty analyses for the ^"''Po a emission rate concentration 



Uncertainty components and propagated uncertainties" 



Relative uncertainties 



(a) LS measurement reproducibility 

(b) Variability in ^**^o a emission rate concentration for any 
ampoule of the SRM 

(c) LS detection efficiency (system and chemical quenching 
(correction) variability) 

(d) Variability in matching LS samples and blanks for background 
subtractions 

(e) Background measurement imprecision 

(f) LS cocktail stability over measurement time (quench variability) 

(g) LS sample volume effect on LS efficiency 

(h) Other LS counting interferences (e.g., wall effects, sample 

heterogeneity, static electricity, etc.) 
(i) Live time determinations for LS counting time intervals 
(j) Variability in LS sample masses 

(k) Gravimetric (mass) determination for any one LS sample 
(1) Extrapolation of LS count rate to zero energy and correction for 

subtraction of 2.3 keV ce response 
(m) Correction of LS response for EC branch decay 
(n) Timing 
(o) Decay corrections for ^***Po from measurement time to common 

reference time 
(p) ^'^^Po impurity correction for embodied uncertainty in ^'***Po/^*^Po 

impurity ratio 
(q) Decay corrections for ^****Po from measurement time to common 

reference time 
Combined standard uncertainty,* « 
Expanded uncertainty' for/: =2, U=ku 



0.012 

0.056 

0.12 (and PE)** 

0.004 

WE** 
0.005 (and PE)** 
0.03 (and PE)** 

WE** 
0.02 
WE** 
0.05 

0.06 

0.06 

negligible (and PE)** 

0.02 (and PE)** 

0.009 (and PE)** 

negligible (and PE)** 
0.17 
0,34 



■ See accompanying text and refer to Sec. 2.L4. 

** The relative uncertainty for this component is wholly (WE), or in part (PE), embodied in the relative 

standard uncertainties of components (a) and (b). 



Component (a), the first in Table 7, is the LS 
measurement imprecision given by the relative 
standard deviation of the mean s^ = 0.012 % for the 
536 independent determinations of /?«. It, of neces- 
sity, contains random variabilities (wholly or in 
part) from many component sources of inaccuracy: 
background subtraction variabilities, spectral inte- 
grations, LS cocktail stability, chemical quenching, 
timing errors, sample mass determinations, and de- 
cay corrections. For comparison, the 5m for 10 repli- 
cate measurements on a single given sample 
typically ranged from 0.05 % to 0.1 %; and 5m for 
all replicate measurements within a given sample 
series was typically ± 0.05 % (see Tables 2 through 

4). 

Component (b) is the uncertainty in mean^a for 
any given ampoule of the ^Po solution standard. It 
was obtained from the standard deviation for the 
between-ampoule dispersion (see Sec. 2.10.1) and 
has V =(7 ampoules — 1) = 6 degrees of freedom. It 
also of necessity contains inherent random variabil- 
ities from other uncertainty components. 



Component (c) is an estimate of the uncertainty 
in the LS detection efficiency, and its variability 
with chemical quenching and over the observed 
system variations (for all of the variables of these 
measurements). The estimate was approximated in 
several ways. The first, call it (cl), was obtained by 
comparing the mean /?a(mean) for all samples and 
measurements within a given sample series (and 
with a given LS spectrometer) to the extrapolated 
intercept value of Ra(ms = 0) (i.e., at zero sample 
mass) obtained from a regression of Ra as a func- 
tion of sample mass nis (as in Figs. 11 through 13). 
Intercept values of Ra obtained from regressions 
with other sample composition variables (e.g., p^q 
or A^aq of Table 1— which also could be used as 
measures of changes in chemical quenching) were 
statistically equivalent to Ra(ms = 0) values. Six 
comparisons of /?a(mean) to RJ^ms = 0), from the N, 
P, and Q sample series obtained with both LS 
counters, had differences ranging from -0.048 % 
to +0.10 %. The mean of the absolute values of 
the sk differences was 0.037 % which may be taken 
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as the component (cl) estimate of the relative stan- 
dard uncertainty. The uncertainty on this uncer- 
tainty (reflected in the mean's sample (n = 6) 
standard deviation was « 90 % on a relative basis. 
Alternatively, the above range of 0.148 % may be 
multiplied by a coefficient of 0.9 (for an unbiased 
estimate of a standard deviation from the range of 
two observations) to obtain a second estimate, call 
it (c2), of M «0.13 %. A third approach (c3) consid- 
ers the difference (range) between the means of 
the most quenched and variable A and B sample 
series (after reduction by the known gravimetric di- 
lution factor) and the mean Aa obtained from the 
n=536 determinations data set (refer to Sec. 
2.10.9 and Table 6). This difference yields, again 
after multiplying by 0.9, a u (c3) estimate of 0.20 %. 
This estimate includes the uncertainty in the gravi- 
metric dilution factor (=-0.05 %). Lastly, for m(c4), 
one can consider the range in^a obtained from the 
difference in the maximum meanv4a (with B series 
samples) to the minimum (Q series samples). 
There are either 9 or 18 mean values to consider in 
the range (see Table 6) depending on if one does 
or does not average the sample series results across 
both instruments. In either case, multiplying the 
range by a coefficient of « 1/3 gives u (c4) —0.12 %. 
On considering these values for m(c1), m(c2), 
u(c3), and u(c4) from the four approaches, a rela- 
tive standard uncertainty of u(c)«0.12 % seems 
reasonable. It also seems reasonable that this un- 
certainty component has in turn a relative uncer- 
tainty of « 70 % (i.e., about 8 parts of 12 for the 
uncertainty on the uncertainty). The uncertainty in 
the LS detection efficiency must also be partially 
embodied in the uncertainty estimators for compo- 
nents (a) and (b) since the latter were obtained 
over a wide range of variables that might have in- 
fluenced the efficiency. 

Component (d) was estimated from the combi- 
nation of comparisons between the quench indicat- 
ing parameters QIP (either H* or tSIE) for the 
samples and matched blanks, and observed func- 
tional relationships between the blank counting 
rates (obtained over narrow time intervals) and the 
QIP, Data for the latter was obtained from a spe- 
cially prepared series of 5 blanks having varying 
sample composition. The dependence of the blank 
counting rate Rh as a function of QIP was generally 
defined in terms of a linear regression having a fit- 
ted slope of dRh/dQIP. The background rate was 
for all intensive purposes independent of QIP, The 
fitted slopes wer€ essentially zero with large inher- 
ent standard deviations: d/?b/dG/P =0.086 ±0.026 



for //'^ and 0.0077 ±0.0086 for tSIE. For a differ- 
ence ^QIP between a sample QIP and that for its 
matched blank, the background counting rate dif- 
ference ARb due to the mismatch is A/?b = 
AQIP(dRJdQIPy The relative standard uncer- 
tainty is then AR\JRh and was typically a few tenths 
of a percent for AQIP differences of ± 1 unit of//* 
or 5 to 10 units of tSIE. The resulting contribution 
of component (d) to the uncertainty inAay requires 
consideration of the sensitivity factor: e.g., M(d) = 
{dAJdRh){^h/Rh)^ The result, based on analyses of 
the P sample series on both LS counters was 
u(d) = 0.0039 % with a 60 % relative uncertainty 
on this uncertainty (i.e., for the 10 determinations 
of ARb/Rh. 

The variability in the background counting rates 
from measurement to measurement (i.e., compo- 
nent (e) uncertainty) was wholly embodied in the 
precision estimators for components (a) and (b). 

The uncertainty due to cocktail stability over the 
measurement time in so far as it affected the chem- 
ical quenching and a detection efficiency (compo- 
nent (f) in Table 7) was estimated by taking the 
variability in the QIP (as given by the standard de- 
viation sqip determined from the multiple measure- 
ments of the QIP for a given sample) and 
multiplying it by the known change dAJdQIP (or 
its equivalent dRJdQIP) in Ra as a function of QIP 
for a given sample series. The estimator Sqip is a 
measure of the variability in QIP over the measure- 
ment time. For example, a linear regression to the 
Ra and //* data for the Q sample series had a fitted 
slope of d/?a/d//* = 0.0057 ±0.0029. For measure- 
ments of the same sample series, the relative stan- 
dard deviation in H* for any sample ranged from 
0.55 % to 0.80 %. The uncertainty in this case is 
then (d/?a/d//*)5/f#, which has corresponding val- 
ues of 0.0015 % at a lower bound and 0.0068 % at 
an upper bound. Analyses of five cases for the N, 
P, and Q series with both spectrometers resulted in 
an overall relative standard uncertainty of 
u(f) = 0.0048 %. Again, as for the component (c) 
and (d) cases, the uncertainty in this uncertainty 
estimate is large, approaching 75 %. The relative 
uncertainty in QIP was Sqip=03 % to 0.9 %. Val- 
ues of dRJdQIP ranged from 0.0001 to 0.006 with 
accompanying relative uncertainties in dRJdQIP 
ranging from 500 % to 50 %. One must further 
note that some part of the random variability in the 
cocktail stability must be included in all of the 
replicates making up the component (a) and (b) 
estimators. 
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The LS detection efficiency varies slightly with 
sample volume and is at an optimal maximum (in 
these sample configurations) for volumes some- 
where between 8 mL and 16 mL. The component 
(g) standard uncertainty estimate was obtained 
from comparing the differences in the mean^a val- 
ues obtained with samples having mtot=10 g to 11 g 
and with that for samples having mtot—14 g to 15 g. 
Part of the variability in LS efficiency due to sam- 
ple volume must also be embodied in components 

(a) and (b). 

Other LS counting interferences (such as those 
arising from discharges of static electricity on the 
LS sample vials or in the LS counter; wall effects 
due to solvent diffusion into plastic vials; sample 
heterogeneity which results in differing efficiency 
in the various heterogeneous phases or compo- 
nents; and random electronic noise) are all spo- 
radic in nature* These effects would increase the 
random variability in the measurement results. 
Considering the wide range of variables over which 
the measurements were obtained (with different 
types of LS vials, over widely varying sample com- 
position, with different instruments and operating 
parameters, and at many different times over more 
than one year), the uncertainty in the LS efficiency 
due to these sporadically occurring interfering ef- 
fects are considered to be wholly incorporated 
within the component (a) and (b) uncertainty esti- 
mators. 

For component (i), the uncertainty in live time 
determinations was estimated by assuming that the 
deadtime corrections that are automatically per- 
formed by the LS counters are appropriately made 
to within ± 10 % (which is a reasonable assump- 
tion given the state of electronics for modern 
ADC-based counting instruments [29]) and on con- 
sideration of the counting rates of the LS samples 
and the known resolving times for the LS instru- 
ments (further assuming that the deadtime inter- 
vals are fixed). The relative uncertainties 
calculated on this basis were 0.01 % and 0.03 % in 
worst cases for the Packard and Beckman instru- 
ments, respectively. Given that roughly half of the 
total number of measurements of /?« were per- 
formed with each instrument, the relative standard 
uncertainty u(i) was taken to be 0.02 %. 

The variability in LS sample mass measurements 
(component (j) uncertainty) is wholly embodied in 
the uncertainty estimators for components (a) and 

(b) in as much as the mean /?„ was based on 27 
independent sample masses (with each mass deter- 
mination having an independent confirmation — see 
Sec. 2.2.3). 



For component (k), the estimated relative stan- 
dard uncertainty in the gravimetric determination 
for any one LS sample is 0.05 % and was treated in 
extenso in Sec. 2.2.3. 

The relative standard uncertainties in the two 
corrections kzE and kEc (for the (1) and (m) compo- 
nents) were respectively addressed in Sec. 2.10.2 
and Appendix A. 

Timing uncertainties were negligible considering 
the long half-lives involved in these measurements. 
For a timing error of 1 min (which is almost incon- 
ceivable), the relative uncertainty in the ^Po decay 
correction is less than 1.3x10"'* %. Also, any tim- 
ing variabilities would certainly have been averaged 
out over the 536 measurements contained in com- 
p)onents (a) and (b). 

For components (o), (p), and (q), the uncertain- 
ties in the corrections for ^Po decay and the ^Po 
impurity (with its respective decay) were derived 
assuming relative uncertainties of u (A^*^Po) = 
4.9 %, u(\^Fo)^0M9 %, and w(/) = 8.1 %. The 
uncertainties in the decay constants A were ob- 
tained from the Evaluated Nuclear Structure Data 
File [30] while that for the ^**^Po to ^**^Po a-emission 
rate impurity ratio / was given in Sec. 2.9. The 
mean Ra was obtained from determinations sepa- 
rated by more than a year (with all determinations 
decay and impurity corrected to an identical refer- 
ence time). Hence, the magnitude of the correc- 
tions (and hence their uncertainties) are quite 
disparate. The 1993 LS measurement results were 
corrected for decay times from 0.7 yr to just over 1 
yr whereas those obtained in 1994 were at most 
only corrected over times of a few weeks (<0.04 
yr). One may therefore consider the ranges of un- 
certainty values over these decay times. For decay 
times ranging from r==^0.04 yr to ==1 yr (see equa- 
tions in Sec. 2.4.1), the propagated relative uncer- 
tainty in the ^'^^Po decay correction ranges from 
0.0013 % to 0.033 %. A mid-value of w(o)===0.02 % 
may thus be used as an estimate of the standard 
uncertainty for this component. Similarly, the cor- 
rection for the ^Po impurity over this same time 
interval ranges from about w(p)— 0.01 % to 
0.008 %. The uncertainty contribution from ^Po 
decay is negligible with an estimated relative 
uncertainty component of w(q)==-8x 10"' % for 
r==0.04 yr to <2xl0-^ % for 7=1 yr. The inher- 
ent random variabilities in these three corrections 
are also of necessity partially embodied in the un- 
certainty estimators for components (a) and (b). 

The treatment used to propagate these compo- 
nents into a combined standard uncertainty and an 
"overall" expanded uncertainty was outlined in 
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Sec. 2.1.4. One should also recall that the data set 
of 536 Ra values that largely comprise these uncer- 
tainties was shown to be normally distributed (Sec. 
2.10.2). 



3. Concluding Notes 

The prepared and calibrated ^^o solution stan- 
dard will be certified and disseminated as NIST 
SRM 4326 with the following specifications: 

Fonn: Solution of nominal 2 mol • L"^ HCl in a 
flame-sealed NIST 5 mL borosilicate glass ampoule 

Solution mass: (5,1597 ±0.0024) g 

Solution density: (1.031 ±0.004) g • mL"^ at 22 °C 

Alpha emission rate concentration: (85.42 ±0.29) 

s-i.g-i 

Known radionuclidic (a emission) impurities: 
(0.106±0.017) S-' • g-' from ^Po 

Reference time: 1200 EST 15 March 1994 

The calibration was based on 4'Tra LS spectrometry 
using two LS measurement systems, with confirma- 
tory measurements by 2ua gas-flow proportional 
counting. The mean a emission rate concentration 
was derived from 536 independent determinations 
of the LS counting rate, using 27 counting sources 
(of widely varying chemical quenching and sample 
composition) that were in turn obtained from 7 
randomly selected sealed ampoules of the ^**^Po so- 
lution standard. The determinations were per- 
formed over a long (=^1 yr) time interval. Photon- 
and alpha-emission impurity analyses were con- 
ducted by 7-ray spectrometry with two intrinsic Ge 
detectors for the former, and by a spectrometry 
with a Si surface-barrier junction detector for the 
latter. 

The 4'Tra LS-measurement calibration of ^**^Po (a 
nearly pure a emitter) would have been routinely 
straightforward otherwise, but was confounded by 
the existence of a low-energy level in ^^Pb (fed by 
the a decay of ^**^Po) that was previously unknown 
to be a long-lived isomeric state. This finding (or 
more accurately, its converse) required the exhaus- 
tive LS spectral investigations that were performed 
as well as the large number of LS measurements 
that were conducted with large variations in chemi- 
cal quenching and sample composition conditions. 



We may ultimately conclude by noting that this 
work is reported here in such depth and exacting 
detail for two reasons. Firstly, because of the com- 
plexity of this calibration; and secondly, and per- 
haps more importantly, to provide a documented 
and archival record of some of the radionuclidic 
metrology procedures routinely practiced by the 
NIST Radioactivity Group. 



4. Appendix A, Correction for LS 
Response to ^^Fo Electron Capture 
Decay 

The LS responses from the counting sources of 
the ^Po standard solution included contributions 
from some portion of radiations arising from the 
^**^Po EC decay mode. Conversion of the previously 
given mean counting rate concentration Ra to a 
^**^Po a-emission rate concentration requires a cor- 
rection for these EC-contributing LS responses. 
The correction can be approximated by consider- 
ation of the known ^°^Po decay and atomic and nu- 
clear structure parameters, and estimation of LS 
detection efficiencies. The following notation will 
be used: 

/a = ^**^Po a decay intensity; number of 
emitted a particles per decay of ^^^''Po 
/ec = ^**^Po EC decay intensity; total num- 
ber of captures per decay of ^Po 
/j = intensity of the p^ radiation compo- 
nent (electrons or photons) arising 
from the EC decay of ^**^Po 
ej = LS detection efficiency for the j^^ 
component 
PK=K-shell EC probability; number of K- 
shell electron captures (holes) per 
EC decay 
6;K==K-shell fluorescence yield; number of 
K X rays emitted per K-shell hole 
PLyPuyP^, .-EC probability for L-shell capture, 
M-shell capture, N-shell ...., respec- 
tively; number of electron captures 
(holes) in the respective shells per 
EC decay 
BK^Bi^yBu,'" = atomic binding energy for the K shell, 
L shell, M shell, .... electrons 
&)L = average L-shell fluorescence yield for 
all three L subshells (including x rays 
following Coster-Kronig subshell re- 
arrangements); total number of L 
X rays emitted per initial L-shell hole 
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6>M = average M-shell fluorescence yield 
for all five M subshells; total number 
of M X rays emitted per initial M- 
shell hole 
a = total internal conversion coefficient 
for a 7 transition; ratio of the number 
of de-excitations which occur by con- 
version (with the emission of an 
atomic electron) to the number of 7- 
ray photons emitted. 

The LS counting rate concentration /?«, when 
normalized to the ^^Po activity concentration, is 
given by L+JJec where /is the fraction of all radia- 
tions from the EC branch that were detected. The 
correction factor A:ec is then 

kEC=IJ(Ia +ffEc) = 1/[1 +f(lEcJIa)l 

The ratio /nc/Za, from the critical evaluation of 
Martin [11], is equal to (0.0048 ±0.0004)/(0.9952± 
0.0004) =0.0048 ±0.0004. The quantity /is a sum- 
mation /=Se// //eg over ally photons and electrons 
arising from the EC decay branch, and where e, 
and Ij/Iec are the LS detection efficiency and rela- 
tive intensity for the p^ component. A lower limit 
on the magnitude of /cec may be established. If all; 
components of the EC branch are detected with 
100 % efficiency (i.e., /=1), then A:ec = 0.9952± 
0.0004 is the lower limit. The upper limit, of 
course, would be A:ec = 1 with/ = (i.e., none of the 
EC radiations are detected). 

Electron capture decay of ^**^Po is a second for- 
bidden unique EC transition (A/'=l/2"->7/2"). 
The result is that K-shell capture is relatively low, 
and that significant EC contributions occur from 
higher atomic shells. This in turn results in the 
emission of lower energy x rays (L, M, ... instead of 
K X rays) and greater Auger electron intensities 
(since ok, (om, ... < (ok). LS detection efficiencies for 
these lower energy x rays and Auger electrons are 
larger that that for K x rays which are usually the 
predominant emissions in high atomic number EC 
transitions. 

The relative intensity of K x rays is IkjJIec=Pk<ok. 
The capture ratio Pl+m+n+JPk was experimentally 
determined by Mandel et al. [13] to be 0.42 ±0.07, 
and which may be compared to a value of 0.39 from 
theoretical estimates of Behrens and Buhring [31]. 
If one adopts a value of 0.4 ±0.1 (and given that 
Pk+Pl+Pm + - = 1), then /'k = 0.70 ±0.015. With 
(OK for Bi taken as =0.970 ±0.015 [32,33], IkJIec = 
0.68 ±0.05. The LS detection efficiency for these 
rather energetic 74.8 keV to 89.8 keV x rays is 



small. Photoelectric absorption is negligible. Their 
interaction in the LS cocktail is primarily by Comp- 
ton scattering which results in Compton electrons 
from zero energy to about 10 keV. As a result, the 
efficiency is estimated to be from near zero to per- 
haps 10 % to 20 %. It will be assumed that 5 % of 
the K X rays are detected, but the upper uncer- 
tainty limit on the detection efficiency will be 
taken as 15 % such that eKx/Kx//Ec = (0.68 ±0.05) 
(0.05l8:i§) = 0.03418:J§4. The relative intensity of K- 
shell Auger electrons is IkJIec =Pk(1 - ^^k) = 
0.021 ±0.002. They have discrete energies ranging 
from about 58 keV to 90 keV given by the atomic 
binding energy differences (Bk-IBu), (^k-2Blii), 
(5k — 2Bliii), (Bk — 2Bmi), ... (Bk — 2Bm\/), (Bk — 2Bsi), 
... etc., and are assumed to be detected with 100 % 
efficiency. Therefore, £kc/kc//ec = 0.021 ±0.002. 

The X rays and Auger electrons for L-shell cap- 
ture can be treated analogously. The relative inten- 
sity of L X rays is IiJIec = {PiJPk)Pkcol* A capture 
ratio of /'iyPK = 0.3±0.1 is taken from the theoreti- 
cal values Pu^JPk = 0.2S and PlJPk = 0.02 [31]. 
The average fluorescence yield for all three L sub- 
shells is a)L = 0.40 ±0.06 which was obtained from 
the mean of six experimental determinations tabu- 
lated by Bambynek et al. [32]. Hubbell et al. [33] in 
a very recent review reported a value of ft>L = 0.382. 
The L X rays have energies of around 11 keV to 15 
keV, and have a range of about 0.6 cm in the LS 
cocktails which is about comparable to the average 
geometric path length in the scintillator volume. 
They are assumed to be detected with an efficiency 
of ei^ = 0.5 ± 0.25, so that eiJIiJIec) = 0.042 ± 0.026. 
The L Auger electrons, having energies of around 
5.4 keV to 16 keV, are assumed to be detected with 
100 % efficiency as for the K Auger electrons. The 
detected fraction of their relative intensity is thus 
eLc(/Lc//Ec) = eUPiJPK)Pyil - a>L) = 0. 126 ± 0.048. 

Higher atomic shell (M + N + ...) EC events have 
considerably lower relative capture probabilities. A 
value of /'M+N+„yPK = 0.08±0.01 is estimated, again 
from theoretical predictions: /*Mi+n/^K = 0.07; PmJ 
/'k = 5.639x10-3; /'m,v//'k = 6.19 x 10-5 [31]. The 
relative intensities for the x rays are relatively low 
and are ignored. Given 6>m = 0.035 ±0.005 [32,33], 
/mx//ec— 0.001 for example. The M shell Auger elec- 
trons range in energy from 0.7 keV to about 3.9 
keV. All N shell Auger electrons have energies of 
less than 1 keV. For convenience (and in the ab- 
sence of additional information with which to make 
better approximations), we can merely assume that 
all of these higher shell Auger electrons are collec- 
tively detected with 50 % efficiency and with a 
somewhat arbitrarily assigned relative uncertainty 
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of ±50 % on the efficiency. With all Auger elec- 
tron yields (l-ci>M), (1 — oin), ..., etc. taken to be 
equal to unity, the detected fraction of the 
relative intensity is therefore assumed to be 

eM + N + ...(/Mc + Nc + .7/Ec) = 0.028 ±0.014. 

Detection of the 896.6 keV 7 ray and its atten- 
dant conversion electron (see Fig. 2) is neglected. 
Firstly, these are in prompt coincidence with the 
radiative and non-radiative EC decays, and there- 
fore would only be detected by the LS spectrome- 
ters if and when the x rays and Auger electrons 
were not detected. Secondly, the LS detection effi- 
ciency for the 7 ray is virtually nil, and the relative 
intensity for the conversion electron, given by the 
transitions total conversion coefficient en, is very 
small. Using 7^ = 0.0047 ±0.0004 and ay = 0.0209 + 
0.0005 from Martin [11], the relative intensity of 
the conversion electron is /cc//ec= 0.020 ±0.002; 
and, as stated, it would only be detected for that 
fraction of the EC events (primarily x rays) that are 
not detected. The undetected fraction of x rays and 
Auger electrons (based on summing all of the de- 
tected relative intensity contributions given above) 
is roughly 75 %. Therefore, the detected relative 
intensity for the conversion electron would be 
<0.015 as a maximum. 

The correction factor may thus be approximated 
by summation of 

^=-(fe)+«-(fe)^-(fc)^-(fe) 

-r eMc + Nc + ...l J j — U.ZM+0.080 

(or/= 0.266 if the maximum contribution from the 
conversion electron is included); which results in 
fee = 0.9988 18:8885. The relative standard uncer- 
tainty on kEc will be taken as the larger upper 
bound, i.e., u = ±0.06 %. 

Note Added in Proof 

Standard Reference Material SRM 4326 has 
been issued with a relative expanded uncertainty of 
0.42 % for the ^**^Po alpha emission rate concentra- 
tion. This is slightly revised from that given here 
(0.34 %). The difference is mainly due to more 
conservative estimates in the livetime measurement 
and LS detection efficiency uncertainty compo- 
nents, and in the treatment of the component 
uncertainties for possible undetected alpha-parti- 
cle-emitting and photon-emitting impurities at 
their limits of detection. Considering the large "un- 



certainty in the uncertainty" associated with these 
estimates, the two values may be viewed as essen- 
tially equivalent. 
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